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Los programas que aparecen en este libro funcionan en los 
ordenadores: 


IBM-PC, XT, AT y compatibles. 


AMSTRAD-464, 664, 6128, 1512. 
SINCLAIR-SPECTRUM 48 K, 128 K, PLUS, PLUS 2. 
MSX-Todos los modelos. 
COMMODORE-CBM 64 y CBM 128. 


INTRODUCCION === 


ASTA la aparición del ordenador, el hombre malgastó mu- 
cho tiempo en realizar operaciones de cálculo y matemá- 
ticas de todo tipo, así como grandes procesos rutinarios. 
El llamado «cerebro electrónico» agilizó todas estas accio- 
nes. Sin embargo, las limitaciones de los ordenadores son 
grandes. 


Queremos dejar constancia de la «torpeza» del ordena- 
dor. A lo largo del libro se subrayará este «defecto» de los 
== microordenadores. Una computadora no resuelve un pro- 
blema, sino que ayuda a su resolución. Hay que aprovechar su mayor cua- 
lidad: la gran velocidad operacional. 


Con la masificación del uso de los ordenadores, más concretamente de 
los personales, el mito de máquina terrible e inabordable parece haber de- 
saparecido. Un micro, mini o gran ordenador no hace más de lo que le 
mandemos; su capacidad de decisión propia es nula, e incluso es de de- 
sear (es una opinión personal) que no llegue nunca a tenerla. 


Lógicamente la aplicación más directa de un microordenador son las 
Matemáticas y la Estadística. Desde que Pascal y Leibniz, con sus ingenios 
mecánicos precursores de las actuales calculadoras, se dieran cuenta de 
la necesidad de descargar este pesado tabajo operacional, en una máquina 
para profundizar más en la Ciencia, continuando por Babbage con sus fa- 
mosas máquinas, diferencial y analítica, hasta hoy día, en que el ordena- 
dor ha invadido nuestras vidas, está clara la utilidad de una herramienta 
capaz de realizar estas monótonas e interminables operaciones, en aplica- 
ciones matemáticas. 


A lo largo de todo este libro abordaremos diferentes temas, en los que 
nuestro ordenador personal puede echarnos una mano. Habrá, natural- 
mente, capítulos en los que la labor del ordenador sea más agradecida, en 
otros su mérito se deberá más al programador. 
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En todo caso se ha tratado de dar un repaso general de materias mate- 
máticas y estadísticas, algunas de ellas poco usuales, procurando valernos 
de nuestro ordenador casero; esperamos que con ello vaya adquiriendo el 
lector el hábito de utilizar su ordenador personal para abordar estas ma- 
terias a partir de ahora. 

Hemos introducido indicaciones, al lector más interesado, de cómo po- 
der ampliar las aplicaciones aquí dadas. Además las matemáticas, quera- 
mos o no queramos, es, quizá, la ciencia más utilizada a diario, pudiendo 
convertir en algo más que una curiosidad el conocimiento de los temas 
tratados. 

Ha sido ardua la tarea de seleccionar, por razones de espacio, sólo al- 
gunos temas entre la amplia gama de posibilidades matemáticas y estadís- 
ticas. La elección ha estado condicionada a presentar un poco de todo, 
abarcando desde los conjuntos hasta las integrales definidas, procurando 
no profundizar demasiado en cada tema tanto matemática como informá- 
ticamente. Este libro está ideado para que pueda ser entendido fácilmente 
por un lector con una base media matemática (e incluso escasa) y con no- 
ciones elementales del microordenador y del BASIC. Los temas son, en ge- 
neral, independientes uno de otro, pudiéndose alterar el orden de lectura. 
De todas formas aconsejaríamos seguir la línea marcada por el libro, para 
su mejor comprensión. 

Finalmente agradecer a amigos y colaboradores su ayuda en la elabo- 
ración de este libro. 


NOTA SOBRE LOS PROGRAMAS 


Todos los programas del libro están realizados con un BASIC 
elemental para que sean válidos para la mayoría de los ordenado- 
res personales. Se emplean siempre instrucciones elementales (IF, 
GOTO, etc.). Están probados en un Amstrad 6128 y en un Spec- 
trum. Además se han revisado manualmente para Commodore. No 
debe dar problemas en cualquier otro micro. 

Sólo se emplean instrucciones gráficas (PLOT, DRAW, etc.) en 
tres programas: estas instrucciones son las que más varían de un 
microordenador a otro, y quizá le exijan hacer alguna pequeña sus- 
titución en su microordenador. De todas formas, se explica el plan- 
teamiento de los diferentes gráficos, facilitando así la labor de 
conversión. 


DIVISIBILIDADERRE 


L uso de las calculadoras ha originado un desuso de las 
operaciones aritméticas básicas realizadas manualmente. 

En este capítulo trataremos una de ellas: la división. 
No vamos a enseñar al ordenador a dividir, esto ya lo sabe, 
sino que intentaremos auxiliarnos con nuestro equipo 
para ciertos cálculos reiterativos en los que intervenga la 
división. Nos referimos a la descomposición en factores 
=== primos, M.C.D., M.C.M. y la aplicación práctica de éstos. 

Los conceptos introducidos son muy sencillos: seguro que usted los co- 
noce. Sin embargo, si no en este momento quizá sí en otro de este libro, 
le resultará divertido repasar estas cifras teóricas para aprovechar al máxi- 
mo todas las experiencias que le proponemos. 

Empecemos por lo más simple, saber si un número es primo. 

Un número es primo, según las Matemáticas, si sólo es divisible por él 
mismo y por la unidad, Sigamos en nuestro programa este planteamiento 
y dividamos el número dado por todos los inferiores a él. Si alguno lo di- 
vide exactamente eso significa que el número que estamos probando no 
es primo. 

Eso es exactamente lo que hacemos en el programa que les ofrecemos 
a continuación: en las líneas 30 a 70 se va incrementando el índice J (lí- 
nea 30) desde 2 (mínimo divisor que probamos hasta N-1 (el mayor) y se 
comprueba si el cociente N/J es entero o no (si N/J es entero, su parte en- 
tera —JNT(N/J)- coincidirá con él y al multiplicar por J en la línea 50 hará 
que R=0). Naturalmente, si N/J es entero, J es divisor de N y N no es nú- 
mero primo. 


Compliquemos un poco este sencillo programa. Calculemos ahora la ta- 
bla de los números primos, desde 1 hasta el número deseado. Observe el 
organigrama 1 


Organigrama 1 


y el programa: 
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Como verá, lo único que hemos hecho es introducir nuestro programi- 
ta dentro de un FOR, para que se repita el proceso tantas veces como in- 
dique el número introducido a modo de límite de la tabla que queremos 
construir. 

Imprimimos el 1 y el 2, que son primos de antemano, para no compli- 
car el funcionamiento de los bucles FOR. 

Pruebe este programa y fíjese en que introduciendo un N muy grande 
el programa se hace muy lento. ¿Cómo ganar velocidad? 

Tanto en este programa como en cualquier otro, hay que evitar opera- 
ciones innecesarias. Si miramos, se observa que no hace falta dividir hasta 
N-1. Sea o no primo, simplemente con repetir el proceso hasta la mitad bas- 
taría. La razón es sencilla: a partir de N/2, es decir, desde la mitad en ade- 
lante, el cociente estaría comprendido entre 2 y 1, nunca poría ser exacto. 
Sólo podría ser exacto en N/2 y, por supuesto, en el extremo superior, N. 
Procedamos a la reforma: 


36 FOR I=2 TO N/2 


Vuelva a probarlo y comprobará la diferencia de tiempo. 
Ya hemos conseguido saber si un número es primo, ahora sigamos el 
mismo razonamiento para descomponer un número en sus factores primos. 


Las instrucciones 40, 50 y 60 ya nos son familiares. En cambio, la con- 
dición del IF ha cambiado. Ahora se interroga si el resto (R) es distinto de 
cero. ¿Para qué? 

Si nosotros directamente incrementamos el índice (J), después de di- 
vidir, en el caso de factor repetido no funcionaría. Lo detectaría la prime- 
ra vez, pero al incrementar, ya no volvería a utilizarlo. Con este error re- 
sultaría que el cociente final (C), que debería ser 1, daría erróneo, o los 
factores resultantes no serían primos. 

Para ello con las instrucciones 66 y 67 volvemos a repetir el proceso 
para un factor primo encontrado, y no incrementará hasta haber agotado 
todas las posibilidades con él. 

¿Cuándo terminará el proceso? 

Matemáticamente sería cuando el último cociente fuese 1, pero esta 
condición no nos bastaría. 

Al despreciar los decimales (INT), podría darse el caso en el que la di- 
visión diese como cociente 1, sin tener en cuenta que, a lo mejor, no fuese 
exacta. Podríamos solucionarlo con: 


Añadiendo la condición R=0, este caso sólo se daría en la última des- 
composición. Pero esto realmente se cumplirá también cuando N<J. Es 
decir, cuando el divisor sea mayor que el dividendo. A partir de este mo- 
mento, las divisiones sucesivas serían inútiles, ya que el cociente (C) val- 
dría siempre 0 y el resto 1. Con lo cual terminaría el programa cuando aca- 
base el FOR. 
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Para el cálculo del M.C.D. y el M.C.M. de dos números nos basaremos 
en dos reglas matemáticas. Podríamos haber seguido el proceso clásico de 
descomponer los dos números en sus factores primos, para su posterior 
estudio. 

Nos encontraríamos con una pega. Los factores habría que ir almace- 
nándolos en algún sitio. La utilización de una tabla sería la solución más 
adecuada, pero estaría limitada a la dimensión (DIM) de la misma. Por tan- 
to, no nos podemos arriesgar a quedarnos cortos ni a desaprovechar 
espacio. 

Utilizaremos dos métodos más rápidos y elaborados, estudiados y de- 
mostrados matemáticamente; debido a ello nos fiaremos de su veracidad. 

Para el M.C.D. el sistema utilizado consiste en dividir el mayor por el 
menor, si el resto de esta división es 0, el menor es el M.C.D. de los dos. 
En caso contrario, entramos en un proceso reiterativo, en el cual repeti- 
mos la operación de la siguiente manera: 

El resto de la división se convierte en el divisor de la próxima y el di- 
visor en el dividendo. Concluyendo, cuando el resto sea nulo, en este mo- 
mento el último divisor es el M.C.D. de los dos números. 

Esto es justamente lo que reflejamos en el programa: 


En las primeras instrucciones comprobamos cuál de los dos números 
es el mayor. A continuación, comienza el proceso reiterativo, que conclu- 
ye cuando R=0 (en la 80). 

En este momento el divisor (M) es el M.C.D. 

Conociendo el M.C.D. conseguir el M.C.M. es muy simple, ya que 


A*B 


M.C.M. = MCD. 
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Incluyamos dos nuevas instrucciones y sustituyamos la 120. 


Ya hemos conseguido saber si un número es primo, descomponer un 
número y conseguir el M.C.M. y M.C.D. de dos números. 

Apliquemos nuestro éxito a una cuestión práctica, que por su elabora- 
ción, no por su dificultad, se hacía pesada. Me estoy refiriendo a la suma 
de quebrados. (Lógicamente la diferencia de quebrados se haría con nu- 
meradores negativos, pero el proceso general es igual.) 


Vamos a calcular la suma de dos quebrados. 
Para ello se seguirá el siguiente proceso: 


1.2 Se pedirán los datos (30 y 40). 

2.2 Se calcula el M.C.M. de los denominadores (50 a la 150). (Esto no 
produce ninguna dificultad, ya que se ha estudiado a lo largo del capítulo.) 

3.2 Siguiendo las reglas matemáticas, obtenemos los numeradores 
(170 y 180). El resultado, lógicamente, será la suma de éstos (190). 


Practique con varias fracciones y compruebe el resultado. ¡Espero que 
le vaya todo bien! 
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Como investigación para el lector podría intentar la suma de varias frac- 
ciones. Observe que no es más que irsumándolas de dos en dos. Como «pis- 
ta» le proponemos una posible solución: sería guardar los numeradores en 
una tabla y los denominadores en otra; sólo haría falta ir incrementando 
el índice e ir calculando la suma por pares. ¡Suerte! 


US 


POLINOMIOS=2 


L estudio de los polinomios y sus posibles raíces se trata- 
rá a lo largo de este capítulo. En todos los programas se 
utiliza la misma filosofía en el tratamiento de los po- 
linomios: 


* Captura de datos: almacenando los coeficientes en 
una tabla de dimensiones N+1 (ya que el término inde- 
pendiente es necesario guardarlo). 

* Tratamiento de un coeficiente: expandiéndolo a to- 
dos los demás gracias a las instrucciones FOR. 

* Impresión del resultado: proceso similar al de la captura, sustituyen- 
do el READ por PRINT. 

Comencemos con lo más elemental: INTRODUCIR LOS COEFICIEN- 
TES Y CALCULAR EL RESULTADO para un valor dado. 


1 


a 


El programa es simple, y en él se aprecian, perfectamente, las tres fa- 
ses explicadas al principio: Captura (10 a 60), tratamiento (65 a la 90) e im- 
presión del resultado (100). 

Superado este paso previo avancemos operando con los polinomios. 
Lo primero será la suma. La regla matemática es sumar los coeficientes 
del mismo grado (200). 


10 REM SUMA DE —POLINOMIOSs POR—FCO. MORALES 


INPUT "INTRODUZCA GRADO DEL POLINOMIO P1 :"3N1 


30 LEFT -Ni=N1+1 


INPUT "INTRODUZCA GRADO DEL POLINOMIO P2 +"5N2 


SÓ LET-N2=N2+1 
60 LET N=5Ni 


IF N2>N1 THEN LET N=N2 


B0—REM-DIMENSTIONAMOS-MATRICES 
90-—DIM P(N)3DIM Q(N)2+DIM S(N) 


140 
130 
160 
170 
180 
190 
z200 


PRINT- "POLINOMIO P1:2" 

FOR J=1 TO N1 
PRINT—"COEFTICTENTE-DE-X*"4 3-14 
INPUT-PCD) 

NEXT-J 

PRINT “POLINOMIO P2:" 

FOR J=1 TO N2 
PRINT-"COEFICIENTE-DE-X*%"5 3-13 
INPUT-QUJ) 

NEXT J 

PRINT:FPRINT "LOS POLINOMIOS QUEDAN ASI:":PRINT 


"POLINOMIO Pis" 


210 
220 
230 
240 
250 
260 
270 
280 
270 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
999 
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FOR J=1 TO Ni 

PRINT PD” 

NEXT-J 

PRINT-"09-=-O"+PRINT 

PRINT "POLINOMIO P2+" 

FOR J=1 TO N2 

PRINT QUID) "xy J=1p9+"“; 

NEXT-3 

PRINT-"0 = O"IPRINT 

PRINT "LA SUMA DE LOS POLINOMIOS ES: " 

FOR J=1 TO MN 

LET-S (EPI 

PRINT-"EL COEFICIENTE DE X“"3J-=15" ES "45(J) 
NEXT J 

PRINT:PRINT “EL POLINOMIO S QUEDARA ASI:":PRINT 
FOR-J=1-TO MN 

PRINT-E (DG JR 

NEXT J 

PRINT “0 = 0" 

END 


Una vez conocida la suma pasaremos a la multiplicación. En ésta, el re- 
sultado se almacenará en una tabla de N+M. Cada elemento de una ma- 
triz multiplica a todos los de la otra, desplazando el resultado tantas posi- 
ciones como las que ocupa. Puede ocurrir que la posición donde vaya a 
dejar el resultado esté ocupada; en tal caso se suman los valores. Este pro- 
ceso viene controlado por los FOR de las líneas 110 y 120. La operación 
y acumulación se da en la 130. 


La división general de polinomios es un proceso más complejo, pero 
puede abordarlo el lector a partir de un caso más sencillo: LA REGLA DE 
RUFFINI. 

Esta se aplica en los casos de divisiones por un polinomio del tipo 
(X +a). El polinomio cociente es de un grado menos; y sus coeficientes se 
calculan (100 a 120) a partir del anterior obtenido (S(J)). El primero es el 
coeficiente de grado N del polinomio (90) y a partir de él se obtendrán los 
demás. 
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Desarrollando el programa (que no presenta ninguna dificultad): 


Nos guiaremos de este programa para el cálculo de las raíces enteras 
de un polinomio. Lo primero que se hace es acotar éstas superiormente, 
es decir, obtener un número que sea superior a todas las raíces (cota su- 
perior). Según Laguerre, un número es cota superior de las raíces de un 
polinomio si al aplicar la regla de Ruffini todos los coeficientes del cocien- 
te y el resto son positivos. ¿Qué haremos? 

Para un valor de A, hay que aplicar la regla de Ruffini y analizar si el 
resultado es positivo (coeficientes y resto). 

Si cumple la condición ya tenemos cota; en caso contrario, incremen- 
taremos hasta conseguirlo. 
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CQHBAY COTA 
0 Lo A A A 


Una vez lograda la cota, si es que existe, intentaremos hallar las raíces. 
El método que utilizan los estudiantes es ir probando número a número 
hasta conseguirlo. Si existe una raíz entera, ésta debe ser divisor del tér- 
mino independiente. Observe el organigrama 2 


Organigrama 2 
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Si se fija tenemos que aplicar la descomposición de un número en sus 
factores primos (TEMA-2), para elegir después los que sean menores que 
la cota. Con éstos probamos uno a uno hasta conseguir hallar las raíces. 
El programa completo quedaría: 


Hacemos un proceso semejante al de búsqueda de los factores primos, 
en el sentido de que si un número es raíz del polinomio puede serlo más 
de una vez. Por ello, cuando detectamos una raíz, «reciclamos» para com- 
probar si vuelve a repetirse. También verá que si probamos con un valor 
A, lo hacemos además con —A. Lógico, ya que —A es divisor del coeficien- 
te independiente. 
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Como aplicación de todo lo dado, sugeriría al lector probar con otras 
reglas para hallar raíces, por ejemplo, las relaciones de Cardano. Biblio- 
grafía de ésta y de otras reglas podrá encontrarse en libros de Cálculo Nu- 
mérico correspondientes a primero de Carrera. Si no se anima, siga ade- 
lante con el libro y no se dé por vencido. 
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BASES DE NUMERACION= 


L estudio de los sistemas de numeración es uno de los te- 
mas más fáciles de comprender y más difíciles de asimi- 
lar. La razón de ello es que desde pequeñitos contamos 
con los dedos de la mano (y de mayores algunos siguen ha- 
ciéndolo) y, por tanto, nos es familiar contar en un siste- 
ma de base 10, pero nos resulta complicado asimilar que 
1+1=10 (en base 2). 

Con la llegada de los ordenadores se empiezan a po- 
===x=x= ner de moda palabras como el sistema binario (base 2), oc- 
tal (base 8) y hexadecimal (base 16). Aquí comienzan nuestros problemas, 
porque el ordenador será muy «hábil» con los 1 y 0, pero nosotros nos ma- 
nejamos mejor con los números del 0 al 9. 

Debido a que nuestro aparato nos ha metido en este embrollo, vamos a 
pasarle el problema y que él lo resuelva (aunque tendremos que ayudarle). 

La base, N, de un sistema de numeración es el número de unidades de 
un orden que se necesitan para anotar una unidad de orden superior. Así, 
en base 10, cuando tenemos tres unidades, por ejemplo, escribimos 3, 
cuando tenemos ocho, escribimos 8 y cuando disponemos de dos veces 
ocho unidades, es decir, dieciséis unidades, escribimos 16 (10 unidades= 1 
decena, y 6 unidades más). * 

De acuerdo con el mismo criterio, en base 6 escribimos el ocho como 
12 (un «grupo de 6»=10 y dos unidades más) y en base 3 el mismo ocho 
lo escribimos como 22 (dos «grupos de 3»=20 y dos unidades más). 

Lo primero es conseguir aclarar el paso de cualquier base a base 10 y 
viceversa. (Vamos a trabajar hasta base 16 por motivos prácticos; si el lec- 
tor quiere ampliar este tope no conlleva gran dificultad. ¡Inténtelo!) 

Empecemos con el paso de base 10 a cualquier base. Matemáticamen- 
te es ir dividiendo el número dado por la base deseada, quedándonos con 
los restos, convirtiendo el cociente en el dividendo de la siguiente división. 

¿Cuándo acabaría? 
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Apoyándonos nuevamente en las Matemáticas, sabemos que un núme- 
ro escrito en cualquier base no puede tener un guarismo mayor o igual 
que la base. Es decir, que si trabajamos en base 4 sólo se puede utilizar el 
OL Zoye3: 

Entonces ya tenemos la condición final; no hará falta seguir dividiendo 
cuando el cociente sea menor que la base. Observe este sencillo programa: 


BACGE NEAR EFA-— MAGA, E 
DAJE Ma FOUR FCO. MORALES 


aS57 ABCDEFGHIJKLM NOPORSTIU 
TR 


ADA AA 


a a A A 
PIDA AS TEN 


Wave 0 ml 
AA e 


Verá una variable, cifra-$, con una serie de valores del O al F. Como di- 
jimos al principio del tema, nos vamos a limitar a 16 posibles bases; por 
tanto, el valor máximo posible será la F. En base 16, como no disponemos 
de símbolos para representar los números por encima del 9, se utilizan las 
letras: A (vale 10), B (valor 11), C, D, E y F (valor 15). Con la función MID$ 
traduciremos los posibles restos en su valor en la base pedida. Lógicamen- 
te hasta la base 10 no haría falta este proceso, pero a partir de ella hay que 
suplir determinados valores por letra. Si el lector quiere ampliar este mar- 
gen de bases sólo tendría que añadir letras a CIFRA$. 

Para el proceso inverso, paso de una base cualquiera a base 10, hay que 
calcular el valor de cada símbolo, según su posición y valor real, acumu- 
lando la suma para el resultado final. Las instrucciones de la 10 a la 30 
son iguales al programa anterior. Aquí el control de fin vendrá dado por 
la cantidad de símbolos que tenga el número, es decir, por la longitud (que 
calculamos mediante la función LEN del BASIC). 

El programa quedaría: 
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50 LET L=LEN(NS$) 

60 -LET-A=0 

70-FOR J=1- TO L 

80 FOR I=1 TO B 

90 TF MIDS$(NS,J,1)<>MID$(C0$,1,1) THEN GOTO 110: REM 
<== EN EL SPECTRUM SUSTITUIR ESTA LINEA POR:+-IF-N 

$(J)<>0$(1) THEN GOTO 110 

100 LET A=A+INT((1I-1)*(B"(L-J))+.5) 

110 NEXT-1 

120 NEXT-J 

130 PRINT "EL NUMERO DESEADO ES: "¡A 

140 END 


Para el paso de una base a otra hay que utilizar como paso intermedio 
la base 10. El programa consistirá en encadenar los dos anteriores, tenien- 
do en cuenta qué base es la inicial y cuál va a ser la receptora del resultado. 


10 REM CAMBIO DE UN NUMERO DE UNA BASE Bi A OTRA BE a 
POR FCO. MORALES 
20 INFUT "INTRODUZCA NUMERO EN CUALQUIER BASE: "¿N 


30 INPUT "INTRODUCIR-SU BASES: ";B1 

40 INPUT “PARA PASARCO-A-BASE: "yB2 

50 LET C$="0123456789ABCDEFGHI JELMNOFQRSTUVWXYZ " 
60 LET L=LEN(N$) 

70 LET A=0 

80 GOSUB 220 

90 FOR J=1 TO L 

100 FOR I=i1 TO Bi 

110 IF —MIDS$(C%1 11 1)<>MIDS$(N4, 3 1) THEN GOTO 130+RE 
M-<=— EN EL SPECTRUM SUSTITUIR ESTA LINEA POR: IF 
C$(1I)<>3NS$ (J) THEN GOTO 130 

120 LET A=A+INT ((I-1)*(B1” (L-J))+.5) 

130 NEXT I 

140 NEXT J 

150 LET-R$="" 

160 LET P=INT (B2x* ((A/B2)-INT(A/B2))+1.5) 

170 LET-R$=MIDS$(0%,P+1)+R$:REM <-— EN EL SPECTRUM 
CAMBIARLO-POR+- R$=CS$(P)+R$ 

180 LET -A=INTCA/B2) 

190 IF A>0 THEN GOTO 160 

200 PRINT "EL NUMERO ";N$;" EN BASE ";B1;" ES IGUAL 
A "R$;" EN BASE ";B2 

210 END 

220 REM COMPROBACION DE DATOS 


22 


Ya le hemos enseñado a nuestro micro a traspasar bases, sigamos edu- 
cándole diciéndole cómo se opera en cualquier base. 

El problema consiste en que tampoco sabemos hacerlo nosotros direc- 
tamente. Para ello nos apoyaremos en lo que sabemos, operar en base 10. 
Ver organigrama 3. 


Organigrama 3 


Este programa lo sofisticaremos un poco, poniendo controles a los da- 
tos. Como indicamos en la introducción, la mayoría de los programas de 
este libro suponen que los datos son correctos. En este caso, vamos a com- 
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izar como ejemplo para 


probar que esto es así y este proceso se podrá util 


los demás programas. 


PORSTIJUVWNYVZ"=] 


OÚDEESGHIJIML] 


UMERO 


A BETARAT uz _ mit - 


| 
| 


DIR DOS NUMEROS (A/B) 


¿O PRINT "4.- DIV 
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conuNTos E 


OR medio de los conjuntos se va a realizar una serie de 
programas, cuya aplicación directa quizá no sea práctica, 
pero valdrán para procesos posteriores. Así, por medio de 
la relación de pertenencia, se tocará el tema de la búsque- 
da, algo muy útil en el tratamiento de datos para seguir 
después con las relaciones de inyección, biyección, etc., 
por medio de las tablas; esto servirá, al lector interesado, 
para el estudio de grafos y circuitos. Terminaremos con 
=== un juego, el MASTER MIND, como aplicación de lo visto 
a lo largo del capítulo. 


¿QUE ES UN CONJUNTO? 


Como todos sabemos, no es más que una colección de elementos. He- 
mos de reseñar que para simplificar los programas, se tratarán conjuntos 
numéricos. Queda a gusto del lector probar con otro tipo de datos. 

Empecemos con la relación de pertenencia, sabiendo que un elemento 
pertenece a un conjunto si está en él. Por tanto, no hay más que buscarlo 
en dicho conjunto. Observe el programa: 


Lógicamente lo primero a efectuar es la entrada al ordenador de los ele- 
mentos pertenecientes al conjunto. 

Después observará que aparte de pedir el dato buscado, se produce un 
proceso reiterativo de lectura-comparación. De esta forma, se lee uno a 
uno cada elemento del conjunto (Tabla), buscando el deseado. Si se en- 
cuentra, ya se sabe que pertenece; en caso contrario, si al recorrer toda la 
tabla no hay éxito, la respuesta sería negativa. 

Este proceso de búsqueda de un elemento es aplicable a los ficheros. 
Como información al lector más avanzado, existen varios métodos. Esto 
se justifica al intentar ganar tiempo en estas búsquedas *. 

Procedamos ahora a analizar el concepto de inclusión. Un conjunto está 
incluido en otro si todos sus elementos están en el otro. Según esto no hay 
más que repetir el proceso de búsqueda para cada uno de los elementos 
del conjunto. 

Fíjese en el programa siguiente: 


SDGEMÉNTAS RE As Mehi 
LEMENTOS DE Az “IN 


SONDNIFa El 
MUDO L£GA 1 


* N. del A.: El estudio de los diversos métodos y la elección para cada tipo de estructura 
es algo muy estudiado por los expertos. Recomendaría, a los interesados en el tema, la lec- 
tura del libro de Algoritmos+Datos=PROGRAMA, en donde encontrará información sobre 
ello. 
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Como verá, no tiene gran complicación este programa. 

Estudiemos ahora las operaciones de la unión y la intersección. Para 
ello se necesita otra tabla que sea el resultado. 

En el caso de la intersección hay que buscar los elementos comunes. 
Aprovechando el programa de la inclusión, se obtendrá fácilmente la 
solución. 

La intersección será repetir la operación de búsqueda de la inclusión, 
con la diferencia ahora de que se irán colocando en una tercera tabla los 
elementos comunes. Esta será la solución: 


La unión la apoyaremos en la intersección. El resultado de la unión se 
almacenará en una tabla cuya dimensión máxima será (N+M). El camino 
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más fácil consistirá en copiar directamente los elementos de uno de los 
conjuntos a la matriz solución. Posteriormente, con el otro se efectúa el 
proceso de búsqueda para cada uno de sus elementos. Si se halla, quiere 
decir que pertenece a la intersección; por tanto, no es necesario repetirlo 
en unión. En caso contrario, se transportan a la tabla unión creada. 


Ya se han hecho las operaciones básicas de la teoría de conjuntos. An- 
tes de seguir con el producto cartesiano y las relaciones de elementos, va- 
mos a crear un programa «sort». 

Un sort es un proceso conducente a la ordenación de datos, ya sea de 
un conjunto, un fichero, etc. Aquí nos vamos a centrar en un conjunto (por 
lógica realmente estamos haciendo la clasificación de una tabla). 

El método más simple, y más largo, es el de la «baraja», así llamado 
por su semejanza al proceso realizado cuando se ordena una baraja de 
naipes. 

Si nosotros cogemos un mazo de cartas para su ordenación, podremos 
ver solamente dos a la vez. Es decir, en un momento determinado se com- 
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para por pares. Fruto de esta comparación, nos quedaremos con la más 
alta (vamos a ordenar de mayor a menor). Una vez elegida la mayor, com- 
paramos ésta con la siguiente, repitiéndose el proceso hasta terminar la 
baraja. 

De esta manera, al final de una pasada completa, está colocada la ma- 
yor en primer lugar. Repitiendo todo el proceso para conseguir la segun- 
da más alta; para ello no empezamos desde la primera, sino ya desde la se- 
gunda. Se reiteraría el proceso N-1 veces. El proceso sería el siguiente: 


MTUNTO ORDENADO QUEDA 
AAA A AA 471 


Desde la instrucción 60 a la 100, se produce el intercambio. Si la carta 
mayor hasta ese momento es mayor o igual que la leída, sigue conservan- 
do su trono. El caso de cartas iguales quedaría de esta forma solucionado, 
ya que sería mayor la que antes se leyera. 

Para obtener el producto cartesiano de dos conjuntos se define una ma- 
triz cuadrada A(N,N), de tal forma que diremos que si dos elementos 1 y J 
están relacionados, el valor A(1,J) es 1; en caso contrario es nulo. 

Este procedimiento de operación es muy útil a la hora de operar con 
la teoría de grafos y en el estudio de circuitos por medio del ordenador. 
Veremos básicamente cómo jugar con ello, pero a los interesados aconse- 
jaría profundizar en el tema. 

El siguiente programa analizará si una relación es de equivalencia. Para 
ello deberá ser: simétrica, reflexiva y transitiva. Estas tres condiciones son 
imprescindibles; en el momento que una de ellas no se cumpliera no ha- 
ría falta seguir. Nuestro micro no sabe nada de conjuntos y como siempre 
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tendremos que proporcionarle el algoritmo correspondiente. Observe el 
programa: 


BO IF CO, 51 THEN PRINT "NO CUMPLE LA PRIPIEDAD REFLEXIVA":GOTO 999 


THEN PRIMT "NO CUMPLE LA TRANSIT 


Entendemos que está bastante claro el proceso. Se ve claramente que 
en el momento que un elemento no cumple una propiedad, el conjunto 
no la cumple y, por tanto, la relación de equivalencia no se produce. Prue- 
be con otro tipo de relaciones. 

Dejando las teorías, juguemos un poco aplicando lo visto anteriormen- 
te. Enseñemos a nuestro ordenador el MASTER MIND o más conocido 
como el juego de los números o de los colores. Consiste en introducir al 
ordenador una combinación de cinco números, y otro jugador tendrá que 
ir probando hasta averiguarla. El ordenador le irá diciendo los muertos y 
los heridos (aciertos en color y posición o aciertos sólo en color sin que 
coincida la posición indicada para dicho color). 


TU ARA E a 
PO PRIMT "INTRODUZCA LA COMBI! 


O 


'O-PRINT- "INTRODUZCA COMEINACION 
(q AFRHLA JAMAL 


Las condiciones adicionales que se han impuesto en el programa son: 
que se juegue con números y que no valga repetir. Si se atreve, pruebe con 
colores (sólo tiene que definir las tablas alfanuméricas) y repitiendo 
colores. 
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COMPLEJOS En 


pesar de lo que su nombre pueda sugerir, realmente los 
complejos no son tan «complejos», sino que las operacio- 
nes necesarias para su manejo predisponen al equívoco. 
En sumas y restas lo más fácil es la forma binómica; en 
cambio, en productos, divisiones y exponenciales lo más 
correcto es la forma polar. El trasiego de operaciones in- 
termedias para trasvasar de una forma a otra origina, mu- 
chas veces, fallos. Todo aquel que haya trabajado con com- 
ZA plejos, sobre todo en su aplicación a circuitos de corrien- 
te alterna, nos comprenderá. 

Se intentará evitar todo este trasiego consiguiendo que nuestro orde- 
nador trate adecuadamente el tema. Nos encontraremos con que nuestro 
equipo comienza a poner pegas. La primera consiste en el uso de las fun- 
ciones trigonométricas; él trabaja en radiones, nosotros más comúnmente 
en grados. Existe un comando en BASIC llamado DEG, que se encarga de 
la misión de pasar de una unidad de medida a la otra. Pero puede ocurrir 
que un determinado modelo no lo tenga (¡no todos son tan listos!). Para 
ello, facilitémosle su tarea, creando unas subrutinas que se encarguen de 
esta tarea. 


==>" REM PASO. DE GRADOS. A RADIANES === 
20 INPUT “INTRODUCIR GRADOS: “¿8 
330 LET R=G*2*3,14159/360 == 
40 PRINT "EN RADIANES ES. ZAR HEEE 


La base de esas subrutinas es una simple regla de tres, teniendo como 
datos fijos 2x1 y 360". 

Ya no tendrá problema con los ángulos; ahora pasemos a operar con 
complejos. Hemos dicho que necesitaremos pasar de una forma binómica 
a polar y viceversa. Programémosle primero para realizar el paso de binó- 
mica a polar, utilizando nuestra rutina de traspaso radiones-grados. 


Se han aplicado las fórmulas de paso de números complejos. El módu- 
lo (30) y la fase (40) se calculan directamente a partir de A y B. El proceso 
inverso es claro, simplemente hay que invertir los papeles según estas dos 
fórmulas: 


e Ni TRIM O ATM 
y r] 1.37] a 


o » A1 


Estamos utilizando la función ATN (arcotangente); si en su modelo no 
existe, debe haber una similar; sustituya ATN por su función. 

Ya podemos hacer todas las transformaciones necesarias. Indicamos al 
principio que para la suma y resta de quebrados la forma binómica era la 
más adecuada. Observe cómo se hace, ya que el complejo suma (resta) no 
es más que la suma (resta) de los respectivos términos de los dos comple- 
jos a sumar (restar). 
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=k 
lO REM SUMA Y DIFERKENC ] 


- = 150 RAS FAM 
0 PRINT “LA SUMA DE LOS DOS COM 


Sencillo, ¿no cree? 
Sigamos con la multiplicación, división y exponenciación: 


Tampoco tiene problemas porque en forma polar es muy fácil realizar 
estas operaciones, utilizando las siguientes fórmulas: 


M1La*M2|$ =M1*M2| 0+B 
Mi1La/M2|B'=M1/M2| a—fB 
(M1iL9)*=M1*L_Xa 


Ya tenemos todo lo necesario para crear nuestro superprograma que 
opere cualquier número complejo, venga en la forma que sea. Para ello bi- 
furcaremos a una determinada dirección dependiendo de la operación a 
realizar, analizando la forma de los operandos, transformándola en la más 
adecuada para operar. 
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3 - INPUT INTRODUZCA PRIMER OPERANDO: aa 


Jr INTRODUZCA SEGUNDO OPERANDO: 


== OPERACION: "*PRINT EEES 


TEA ERA EZRA 


"2.— RESTA" ===== 
3. MULTIPLICACIÓN" 
A DIVISION" = AS 


9 = 
--400-ENPUE-D- -——--- AA == = — 
= 110 IF O=1 THEN GOTO 200. : => 


117-1F G=2 THEN GOTO 300 


== 420 END 

200 1F A$="P" THEN GOSUB 600 

210 LET AS=A1+A2 

= 220 LET BS=B1+B2 

230 PRINT "LA SUMA EN FORMA BINOMICA ES: 
240 GOTO 999 

300 1F As="P" THEN GOSUB 600 


320 LET BS=B1-BZ 
330 PRINT "LA RESTA EN FORMA BINOMICA ES! 
340 GOTO 999 
400 IF A$="B" THEN GOSUB 700 
410 LET ASSAIXAZ 
420 LET BS=B1*B2 


—"AS,BS 


; "yAS/BS 


430 PRINT “LA MULTIPLICACIÓN EN FORMA POLAR ES1 " 


Es 
440 GOTO 999 
-500 IF A$="B" THEN GOSUB 700 
510 LET AS=A1/A2 : 
520 LET BS=B1/B2 
== 530 PRINT “LA DIVISION EN FORMA POLAR ES: 
540 GOTO 999 
600 LET Mi=A1:LET M2=42 
10 LET Al=M1x*COS(B1)1LET Bi=M1*SIN(B1) 
620 LET A2=M2*C0S(B2):LET B2=M2*SIN(B2) 
630 RETURN 
-700 LET Xi=A1:LET X2=A2 
710 LET A1=SQR(X1"2+B1”2) 
720 LET A2=SQR(X2"2+B2"2) 
730 LET B1=ATN(X1/B1) 
740 LET B2=ATN(X2/B2) 
750 RETURN 
999 END 


"AS,BS 


Hemos dejado para lo último la radicación de un complejo. El motivo 
es que aparece un fenómeno curioso en esta radicación. Con complejos, 
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la raíz N da lugar a N raíces, separadas 360%/N una de otra. El programa 
calcula el valor de estas raíces hasta un N; fíjese cómo el módulo y el án- 
gulo de fase disminuyen, ya que 


VMLa) = YMLa/N. 


"A COMPLEJO EN FO 
LAA ATTE EI TINTA 


Este programa gráficamente es bastante interesante. Como las ins- 
trucciones gráficas de un micro a otro varían, vamos a indicarle cómo 
hacerlo: pruébelo. Según aumenta N, las circunferencias formadas por 
las raíces van siendo más pequeñas y los radios que unen el centro con 
las raíces se juntan. Dibuje una circunferencia por cada valor del radio 
(R=/módulo), y coloque los puntos de las raíces; tenga en cuenta que 
tendrán que ir girando 360%/N. Pruébelo, le gustará como queda. 
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POLIGONOS 6 


STAMOS consiguiendo progresos en nuestro micro; inten- 
temos ahora introducirle las fórmulas de geometría para 
descargar nuestra memoria, ya que la suya funciona me- 
jor. No por capacidad; la capacidad del hombre no es com- 
parable a ninguno de los grandes monstruos de ordena- 
dores, sino por estar más «despejada». 

Juguemos primero con los pequeños polígonos, cua- 
drados, trapecios, rombos y rectángulos (los triángulos los 
trataremos en el siguiente capítulo). 

Hagamos que el ordenador nos vaya interrogando hasta tener los datos 
necesarios para el cálculo del área. Esta se calculará aplicando las super- 
conocidas fórmulas geométricas respectivas; si no las recuerda en este mo- 
mento, las encontrará en las líneas 110, 210, 310 y 410. 


A 
¿21 THEM INPUT “INTRODUZCA LA 


SA THTAGANA! 
INTRODUZCA DIAGONAl 
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Vamos a «jugar» ahora con el círculo. En él se dan conceptos (segmen- 
to, corona, etc.) no tan habituales y utiliza fórmulas en el programa para 
que nuestro ordenador aclare la situación. Vamos a emplear un «MENU»: 
que es un buen método para elegir opciones. También se utilizan las ins- 
trucciones ON GOSUB y ON GOTO para bifurcar según la opción teclea- 
da. El ordenador nos irá pidiendo los datos necesarios a medida que los 
necesite; al igual que en el anterior programa, encontrará las fórmulas apli- 
cadas en el propio programa. 


Como se observará, se ha utilizado la técnica de ir preguntando al usua- 
rio del programa por los datos necesarios, en vez de exigirle que los intro- 
duzca. Esta técnica es muy útil cuando van a utilizar el programa personas 
que no conocen el ordenador ni el proceso seguido en su programación. 

A continuación vamos a ver cómo calcular el área de un polígono to- 
mando como dato el radio de la circunferencia circunscrita. Si mantene- 
mos fijo éste y aumentamos el número de lados, resulta que el área del po- 
lígono se acerca al de la circunferencia. Si incrementamos mucho el nú- 
mero de lados llegará un momento en que el polígono y la circunferencia 
se confundan. Veamos esto con el siguiente programa. Imprimimos cada 
N veces para no hacer muy grande la solución. 
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Fíjese en el resultado obtenido por nuestro programa. 


Resultado para 

n=5. 

Este vértice es fijo 
a partir de él se 
acen los lados. 
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Esta aproximación se nota más gráficamente. Como sucedía con los 
complejos, el procedimiento es dibujar la circunferencia, y colocar los vér- 
tices girando con el ángulo 360%/N; después unir dos a dos éstos. El pro- 
grama resultante lo encontrará en el programa siguiente: 
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3O INPUT "INTRODUZCA VALOR INICIAL: "zT1 

40 INPUT "INTRODUZCA VALOR TOPE: “;T2 

50 SCREEN 2:REM *SCREEN 2*<-— VALIDO PARA M.S.X.>» 

PARA EL 1.B.M. CAMBIARLO POR *SCREEN 1?, Y QUITARL 
O PARA SPECTRUM Y AMSTRAD 

60 FOR J=0 TO 2*3.14159 STEP .07 

70 PSET (70+70*C0S(J),70+70*SIN(J))3REM *PSET* <-- 

VALIDO PARA M.S.X. e I.B.M., PARA LOS DEMAS CAMBIA 
R EL *PSET*? POR *PLOT?* 

80 NEXT J 

90 FOR I=T1 TO T2 

100 LET A=70+70*C0S(0):LET B=70+70*SIN(0) 

110 PSET (AsB):REM *PSET*<-- VALIDO PARA M.S.X. e 

1.B.M., PARA LOS DEMAS CAMBIAR EL *PSET* POR *PLOT 
, 


120 FOR J=1 TO I 

130 LET X=70+70*C0S (6. 2832/1+*J) 

140 LET Y=70+70*SIN(6.2832/1+J) 

150 LINE -STEP (X-A,Y-B):REM *LINE -STEP*x<-- VALID 

O PARA M.S.X e I.B.M., EN EL AMSTRAD SUSTITUIR POR 
"DRAWR? Y EN EL SPECTRUM POR *DRAW?” 

160 LET A=X 

170 LET B=Y 

180 NEXT J 

190 NEXT I 

200 GOTO 200 
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TRIANGULOS L 


N el capítulo anterior vimos la manera de calcular las 
áreas de los distintos polígonos. No abordamos en él los 
triángulos, debido al tratamiento especial que le vamos a 
dar a este tema. Esta figura es la más estudiada por los ma- 
temáticos, y base de los cálculos realizados en las demás. 
Trataremos las funciones trigonométricas; éstas represen- 
tan relaciones existentes entre los lados y ángulos de un 
triángulo. 

=== Empecemos creándonos nuestra propia subrutina para 
el cálculo del seno. Nos basaremos en el desarrollo de Taylor: 

Xx? x? xQn+1) 


=Ya RE A AAA 
a aia REE y E 2n+11 


En esta expresión obtendremos mejor o peor aproximación dependien- 
do de la cantidad de términos que tomemos; cuanto mayor sea ésta, más 
exacto será nuestro cálculo. Para nuestro proceso necesitamos primero 
calcular el factorial de un número; lo veremos levemente, ya que en el ca- 
pítulo 12 lo trataremos con profundidad. 


Lo interesante de este programa, a partir del cálculo del seno, es: 


1.2 Conocer el factorial de un número. 
2.2 El manejo de un proceso reiterativo controlando el final. 
3.2 La creación de subrutinas que faciliten la resolución de un pro- 
blema complejo. 
A continuación aparecen el programa y el organigrama 4. 


NV PAL AMA AENA a a 
ALCUCLO DEC VA j a 
an ” 


30 SaSUa 00 DI NDA CAC (ZRN+1))/ FA) 
70 1F ABS(SN-AN)<SE-05 THEN PRINT "EL VALOR BUSCAD. 


Organigrama 4 


DL 


Como observará, el programa se compone de un proceso reiterativo ba- 
sado en Seno=Seno+(—1)""! * X?2"+1/FACTORIAL, siendo el segundo 
seno el acumulado hasta el momento. El error vendrá dado por la diferen- 
cia entre el hallado y el calculado anteriormente. Si esta diferencia es mí- 
nima (<0,005), ya no es necesario seguir. La subrutina factorial (300-360) 
se estudiará más adelante. 

Comenzaremos con el triángulo, rectángulo dada su simplicidad, gra- 
cias al ángulo recto que tiene. Existen cinco posibles incógnitas, dos án- 
gulos y los tres lados; necesitamos dos para conocer las demás, una de és- 
tas debe ser un lado. La razón es lógica, conocido los tres ángulos sólo po- 
demos saber la proporción entre los tres lados; haría falta un lado para 
calcular los otros dos. Haremos que nuestro micro nos vaya pidiendo los 
datos iniciales, él después hará los cálculos. Existe una primera parte (10 
a la 96) en la que se interrogará sobre la clase de datos disponibles, bifur- 
cando, según éstos, a la dirección de programa adecuada. En los diferen- 
tes puntos a los que se ha bifurcado se pedirán los datos a aplicar en las 
fórmulas. Como en todos los programas de geometría, los procesos de 
cálculo no tienen dificultad alguna, pues se limitan a incluir en la línea de 
programa correspondiente la fórmula matemática que corresponda apli- 
car. : 


lO REM FALFINMN E AREZ 
LY TMREPT CALCULO DE AREA FARA 1 


200 INPUT "INTRODUZCA | 
214 —— AMIA 1 


No hemos tenido inconveniente en abusar de IF y GOTO, para simpli- 
ficar el razonamiento. Bifurca por un MENU y por un ON GOSUB o un 
ON GOTO. Acostúmbrese a razonamientos simples para mejorarlos con las 
herramientas que su equipo le permita. Tenga en cuenta que los IF y GOTO 
existen en todos. 

Como en el tema 15 (Espacio vectorial CF) veremos estos conceptos, 
ahora investigaremos sólo las alturas para el cálculo de las áreas. Añadien- 
do al programa anterior estas instrucciones, obtendremos el valor del área. 
Ahora bien, como el que trabaja es nuestro micro, aprovecharemos para 
que nos compare los tres posibles del área; según se calcule con una u 
otra altura, lógicamente deben coincidir. 


Si introduce cualquier valor, comprobará cómo coinciden las áreas; en 
caso contrario, es que hay algún error: verifique sus procesos. 
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RESOLUCIÓN 
DE ECUACIONES (1) 


RATAREMOS de conseguir que nuestro ordenador sea ca- 
paz de resolver ecuaciones de primer grado de una y dos 
incógnitas. Lo iniciaremos con algo fácil: sacar una tabla 
de valores según su ecuación. 


Si ejecutamos el programa anterior en nuestro microordenador, obten- 
dremos los siguientes resultados: 


———————— 


ul 
0] 


3 2 
4 4 
S 6 
6 8 
7 10 
8 E 
s 14 
10 16 


Supongamos que queremos saber para qué valor de X sucede que Y=0; 
es decir, cuál es la solución de la ecuación. Lo más fácil será resolverla des- 
pejando la X: 


B 
cc 


No tendría ningún mérito introducir esta fórmula en nuestro or- 
denador. 


Lo primero, porque hay que buscar una solución válida para otros ti- 
pos de ecuaciones, no simplemente para una recta. Y lo segundo, porque 
no tendría sentido alguno hacer trabajar nuestro ordenador para tal 
menudencia. 


Si el coeficiente de la X no es igual a O, la recta tiene que cortar en al- 
gún punto al eje X. Este punto de corte es la solución de nuestra ecuación. 
En este punto la recta queda dividida en dos partes: una positiva y otra ne- 
gativa. ¿Cómo solucionarlo? 


Para ello nos basaremos en un método matemático que consiste en to- 
mar dos valores de la recta, uno positivo y otro negativo. Esto querrá de- 
cir que la solución se encuentra entre estos dos. Dividimos el intervalo for- 
mado por estos dos puntos en dos partes iguales. Calculamos el valor del 
punto medio; si es negativo, tomaríamos el intervalo formado por el ex- 
tremo positivo y el punto medio. Si, por el contrario, fuese negativo, el pro- 
ceso sería inverso. De esta forma, nuestra solución se encontraría en este 
intervalo, con lo cual hemos reducido a la mitad el número de puntos. Si 
repetimos este proceso llegará un momento que el intervalo quedará re- 
ducido a un punto, ésta es la solución. 
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Si tuviéramos que hacerlo a mano, nos volveríamos viejos operando. 
Utilicemos a nuestro ordenador como un trabajador infatigable. El orga- 
nigrama resultante será: 


Organigrama 5 
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Traspasando el organigrama 5 a programa: 


y 


Quizá no vea práctico este programa para una ecuación tan sencilla, 
pero si lo aplica a una función más complicada podrá obtener los puntos 
de corte con el eje X. También es interesante la selección de un intervalo 
u otro, según el valor del punto medio; esto le ayudará a la hora de elec- 
ción de un camino en un proceso reiterativo. 

Notará que preguntamos el intervalo inicial; lógicamente, si estos dos 
datos iniciales dan valores de Y del mismo signo, no nos sirven; habrá que 
introducir otros: nuestro micro no puede saberlo todo. 

Pasemos al caso de dos ecuaciones con dos incógnitas. Lo resolvere- 
mos de dos maneras. La primera será programar al ordenador para que 
sea capaz de utilizar los métodos tradicionales: reducción, sustitución e 
igualación. Lo prepararemos para reducción. Para crear la lista de instruc- 
ciones seguiremos los pasos matemáticos del método. Fíjese en el organi- 
grama 6 y programa correspondiente: 


Organigrama 6 


Para el segundo método nos basaremos en el razonamiento hecho para 
una incógnita. Indicamos que la solución del sistema consistía en dividir 
la recta en dos partes de distinto signo. Aquí la comparación es entre dos 
rectas. Antes del punto de corte, los valores de una recta son mayores que 
la otra; a partir de la solución, se invierten los papeles. Aprovechando la 
deducción del método para una incógnita, comprobaremos no el valor de 
Y, sino el valor de la diferencia de los valores de Y. El porqué se debe a 
lo dicho anteriormente: si restamos los valores de Y antes del corte, las di- 
ferencias serán todas negativas o positivas, y a partir de éste, tomarán los 
valores contrarios; luego el razonamiento es válido; tendremos que ir par- 
tiendo los intervalos hasta conseguir dar con la solución. Para ver mejor 
el proceso añadiremos una instrucción PRINT que nos vaya informando 
de los valores con los que va trabajando el ordenador. 

Nuestro programa definitivo: 


== REM CALCULO: DEL PUNTO DE CORTE DE DOS. RECTAS. 
15 DEF. FNA(X)=2x*X-4- = == 
16 DEF FNB(X)=X+1 ==> 
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(F] (XIII ENI 


— YO 4 a a 
PAT ATI NALI FEINA ISA 


Probémosle con un intervalo pequeño y nos devolvería los siguientes 
datos para dos rectas fáciles: 


Y=2X-4 
Y=X+1 


En el siguiente programa vienen los resultados que nos devuelve el 
ordenador: 


Esperemos que con estos procesos descritos haya quedado claro cómo 
conseguir resolver ecuaciones de primer grado. 

En sí lo importante no es el resultado, sino el desarrollo empleado para 
conseguirlo. En temas posteriores resolveremos un sistema de ecuaciones 
con N incógnitas. 


6 


a 


RESOLUCION 
DE ECUACIONES (II) 


RACIAS a nosotros, nuestro ordenador va mejorando y ya 
sabe resolver ecuaciones de primer grado. Instruyámosle 
más para que sea capaz de resolver ecuaciones de segun- 
do grado. 

Para la resolución de estas ecuaciones no vamos a se- 
guir el mismo método que para las de primer grado. Nos 
fundamentaremos en la fórmula general de resolución de 
ecuaciones de segundo grado. Esta es: 


X=-—B+yB* 4AC/2A 


Empezaremos por preparar una subrutina que, dados los valores a, b, 
c, nos devuelva las dos soluciones X1 y X2. Numeramos las instrucciones 
a partir de la 500, para poder llamar a esta subrutina en posteriores pro- 
gramas. La primera instrucción (la 10) nos sirve para poder emplear esta 
subrutina como programa independiente. La línea 600 (END) deberá ser 
sustituida por un Return, a la hora de introducir la subrutina en un 
programa: 


“EZÓ PRINT 11Ac 
MD TFT] LAS 


Como ve, no tiene dificultad alguna. Pruébela con distintos datos y com- 
pruebe las soluciones. 
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Al calcular el discriminante (b?—4ac) puede ocurrir que sea menor que 
O. Si se deja así, al ser negativo, nuestro ordenador se quejaría en la fun- 
ción SQR (línea 70), dando error e interrumpiendo el proceso. Es dema- 
siado torpe para mandar un mensaje y continuar la ejecución. Necesita 
una ayuda por nuestra parte. Como siempre, el hombre es superior a la má- 
quina. Tendríamos que retocar la estructura del programa. Analicémoslo 
con el organigrama 7, donde se examinan todos los posibles casos. 


Organigrama 7 


El planteamiento estudiado en el organigrama es el correcto; están to- 
madas en cuenta todas las posibilidades (incluida la solución compleja). 
En este último caso se producirían dos soluciones conjugadas, que ven- 
drán dadas en las líneas 510 y 520. En los otros dos casos utilizamos las 
fórmulas expresadas en el organigrama. 
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O INPUT "INTRODUZCA Ay,By/! 
Se ¿MIRNA Hy y 


Pasemos a ver un caso particular de las ecuaciones de segundo grado. 
Son las bicuadradas. Estas, realmente, son de cuarto grado, pero por su si- 
militud y tratamiento se solucionan basándose en las de segundo grado. El 
razonamiento es fácil: manejarlas como si fueran ecuaciones de segundo 
grado y, posteriormente, hallar la raíz cuadrada de cada raíz. 

No existe dificultad alguna, exceptuando el caso en el cual las raíces 
sean complejas. Observará que el programa emite dos soluciones y un 
mensaje (70 a 90); estas soluciones son las calculadas al aplicar el proceso 
de una ecuación de segundo grado; posteriormente, habría que aplicar la 
radicación de complejos (capítulo 5), obteniendo las verdaderas solucio- 
nes de la bicuadrada. 


JT "Yi=""*EMRA 
IN ms : 


B=-SORCBOS-24xXA6%0) 
AA dd 


Esperemos que no haya tenido ninguna duda. 

Existe otro tipo de ecuaciones: las trigonométricas. Son aquéllas en que 
la incógnita no es X, sino una función trigonométrica de X. El proceso es 
similar en las bicuadradas. Se trata la función como si fuera de segundo 
grado. Con las soluciones se calcula el arco que corresponda. Por ejemplo 
(si es Tg(X)), sería ATM(X1) y ATN(X2). 

¡Intente hacerlo solo! 
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MATRICES 1 


N el tratamiento de las matrices está especialmente indi- 
cado que intervenga el ordenador. La explicación es muy 
lógica; los procesos son simples, pero enormemente te- 
diosos y en este tipo de procesos es donde se obtienen las 
mayores ventajas de un ordenador, por su rapidez de cálcu- 
lo. Los procesos matriciales se centran siempre en un ele- 
mento determinado; posteriormente se aplica a cada uno 
de los distintos elementos de la matriz. Por tanto, todo se 
reduce a programar al ordenador, resolviendo el proble- 
ma concreto de uno de los componentes de la matriz para después, por me- 
dio de los FOR'S, extenderlo a los demás. 

Aplicando lo expuesto, calcularemos la transpuesta de una matriz. Sim- 


E 


FA SFM AMAIA AIM FIAR 
¿DVD DIM AMAN ADM TAM N) 
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160 PRINT T(Jk) 
== 


180 NEXT J EZ = === == AR RA —= 
ÁS a : == == 


En este programa se aprecian las tres fases: recogida de datos (40-90), 
aplicación a un elemento (80) e impresión (120 a 170). No tiene ninguna 
dificultad, ya que sólo hay que igualar el elemento de la matriz a otro, cu- 
yos índices estén intercambiados. 

Continuemos operando con estas matrices. La suma y diferencia de dos 
matrices se obtiene sumando, o restando, un elemento con su correspon- 
diente en la otra. Lógicamente deben tener las mismas dimensiones. 

El programa resultante quedaría: 


10 
20 


REM CALCULO DE LA SUMA 
INPUT "INTRODUZCA LA DIMENSION DE LAS MATRICES: 


"5N,M 


30 


DIM A(N,M)>:DIM B(N,M) 
PRINT "ELEMENTOS DE LA MATRIZ NUMERO 1:" 
FOR J=1 TO N 
ERIMNESPITEASSST 
FOR K=1 TO M 
INPUT A(Jy,K) 
NEXT K 
NEXT J 
PRINT "ELEMENTOS DE LA MATRIZ NUMERO 2: " 
FOR J=1 TO N 


140—PRINT-"FTIEA"33 
150 FOR K=1 TO M 


INPUT B(J,K) 

NEXT K 

NEXT J 

PRINT “LA SUMA ES:" 
FOR J=1 TON 

PRINT "FILA "yJ 

FOR K=1 TO M 

PRINT A(J,K)+B(J,K) 


240 NEXT K 


NEXT J 


Queda para el lector insertar instrucciones que controlen posibles erro- 
res en los datos, tanto de dimensión como de elementos. 
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Como prueba de lo indicado al principio del tema, sólo hay que estu- 
diar el problema para un elemento (instrucción 230), aplicando los FOR'S 
para operar con todos los términos. 

La multiplicación de matrices constituye un caso interesante del mane- 
jo de tablas. Aquí cada elemento resultante es suma de los productos de 
los términos de la fila 1 con los de la columna 1 de la otra matriz. Pasan- 
do al programa: 


10 REM MULTIPLICACION DE MATRICES 

20 INPUT "INTRODUZCA DIMENSION DE MATRIZ 1: ";F1,C1 

30 INPUT "INTRODUZCA DIMENSION DE MATRIZ 2: ";F2,C 2 

40 IF C1<>F2 THEN PRINT "NO ES POSIBLE LA MULTIPLI 
-CACION":GOTO 999 = 

45 DIM A(F1,C1):DIM B(F2,C2):DIM P(F1,C2) 

50 PRINT "ELEMENTOS DE LA MATRIZ NUMERO par 

60 FOR J=1 TO Fi 

70 PRINT "FILA ";J. 

80 FOR K=1 TO C1 

90 INPUT A(J,K) 

100 NEXT K : 

110 NEXT J : y 

120 PRINT "ELEMENTOS DE LA MATRIZ NUMERO. 2. 

EA EEE == 

140 PRINT "FILA "Y 
- 150 FOR K=1 TO C2 
- 160 INPUT B(J4Kk) 
170 NEXT K 222 
180 NEXT J——-- 

190 FOR J=1 TO Fi SS 

200 FOR K=1 TO C2 == 

210 FOR L=1 TO C1 - 

220 LET P(J,K)=P(J,K)+A(J,L)*B(L,K) 

230 NEXT L 

A AAA 

250 NEXT J == 

260 PRINT "LA SOLUCION ES:* 

270 FOR J=1 TO F1 = = 

280 PRiNT "FILA "5J - 

290 FOR K=1 TO C2 
300 PRINT P(J,k) === 
== AAA 


Lo más interesante de este programa es la anidación de bucles; es de- 
cir, el incluir un bucle dentro de otro y en su interior otro, etc. (190-250). 
Así se controlan las filas y columnas de las matrices pudiendo aplicar la 
fórmula de la multiplicación en la 220. 
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En este programa es necesario el control de la instrucción 40, debido 
a la condición imprescindible de C1=F?. Sin este requisito el tercer FOR 
(línea 160) no controlaría todos los elementos, originando errores no 
deseados. 


Para el cálculo de la matriz inversa es necesario saber el valor del de- 
terminante. Se plantea el programa basándose en las propiedades de los de- 
terminantes. Si tenemos una matriz A de orden N tal que: 


a11 212 --- A1n 
A= [271:A32 00: Az 
An1 An2 +++: Ann 
Debemos aprovechar las herramientas matemáticas para lograr trian- 
gularizar esta matriz de forma que quede: 


paa 0 a, 
ll 
0 0 a 


¿Cómo conseguir que nuestro ordenador lo haga? 

Centrémonos en un elemento perteneciente al triángulo inferior de la 
matriz. Este elemento debe convertirse en 0; para ello aplicaremos la si- 
guiente propiedad: 

Si una matriz B se forma sumando a una fila cualquiera de A otra fila 
de A multiplicada por un número real: 


BJ] =[A] 


Sea Aij (¡>j) perteneciente, lógicamente, a la fila i-columna j. Si a Aij le 
restamos 


As 
Ma 


se nos anula. Pero esto mismo, según la propiedad expuesta, ha de apli- 
carse a toda la fila i. Reiterando el proceso a toda una columna se nos anu- 
larán todos los elementos Aij (i>j). Planteando el organigrama 8: 


Observe cómo se controla si el determinante es 0. Este caso se produ- 
cirá cuando los elementos de una columna, a partir de la diagonal inclui- 
da, sean nulos. 


También hay que tener en cuenta el posible inconveniente de que sea 
nulo un elemento de la diagonal; en este caso se aplica otra propiedad. 
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Organigrama 8 


Si una matriz B se forma a partir de una matriz A mediante el inter- 
cambio de las filas cualesquiera: 


¡BS DA 
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El programa resultante, teniendo en cuenta todas estas consideracio- 


nes, es el siguiente: 


Se observan varias partes: 


1.2 De la 60 a la 100 se introducen los datos; este proceso es similar 


a los demás programas del capítulo. 


L Zo 


2. De la 130 a la 170 se comprueba si es nulo el determinante. 
3.2 A partir de la 180 se produce el proceso de triangulación. En la 
220 se halla el valor de S(C), para después aplicarlo a cada fila utilizando 
como coeficiente a multiplicar el conseguido en la 240. Este se extenderá 
a los demás elementos de la fila en el FOR comprendido desde 290 a la 320. 


Para el cálculo de la inversa el proceso es idéntico; simplemente aña- 
diremos unas líneas de programa, de tal forma que los procesos efectua- 
dos a la matriz se reflejen en la matriz identidad. Al final, la matriz resul- 
tante de estas operaciones es la inversa. 


y -—INI = DI RR BNIMENCSIM 
A . AA « — SA y 


A A A A A 1 o O a 
¡o ¡92 =LES L1iCsII=Z1 (FILE 1(F9))=L] 


13 
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SISTEMAS DE ECUACIONES 
DE N INCOGNITAS 


RACIAS al desarrollo de los procesos matemáticos matri- 
ciales descritos en el capítulo anterior, se puede resolver 
con relativa simplicidad el problema de los sistemas de N 
ecuaciones con N incógnitas. Dado un sistema represen- 
tado por: 


aX, +a¡,X7+a¡,X3+ ... +ajX,=C, 
aX, +a,X,+a2X,+ ... Fa X, =C, 
aX, +azsX,+ayX,+.. +a,yX,=C; 


211 A] Tas X Ha, Ag. Fan, Cr, 


Conseguir los valores X1, X2, ..., Xn que satisfagan la igualdad es tarea 
para nuestro ordenador. Para ello, se tratará de construir programas a ima- 
gen y semejanza de los métodos matemáticos existentes. 


Hay que tener cuidado en controlar con precisión los bucles necesa- 
rios para el desarrollo de los métodos. Es importante no confundir los ín- 
dices de las FOR, sabiendo en todo momento qué clase de variable repre- 
senta cada uno, ya sea el número de la columna o el de la fila. 


Antes de empezar convendría tener en cuenta una serie de propieda- 
des utilizadas en los razonamientos. Estas son: 
e Un sistema no varía si se intercambian dos ecuaciones entre sí. 


e Se puede multiplicar todos los términos de una ecuación por un nú- 
mero real distinto de O. 


e Estudiaremos dos métodos conocidos: el de Gauss y el de Gauss-Jor- 
dan. Empecemos con el segundo de ellos. 


TS 


METODO DE GAUSS-JORDAN 


Este se basa en conseguir un sistema de la forma: 


a,¡X¡+a19X,+29/X3+ ... aX, =C, 
aX, +ta,X3+ ... Fa, X,=C, 


e O. AA E O a1=Ca-1 


Si se tuviera un sistema de este tipo, la solución sería trivial. Partiendo 
de la última ecuación se despeja X,.: 


Llevando este valor a la penúltima se calcula X,,_,, y repitiendo este pro- 
ceso de abajo a arriba, se logrará hallar los restantes valores deseados. 
Pero el conjunto del que partimos no es este sistema triangular, sino 


el siguiente: 


aX +a¡,¿X¿+a¡¿X3+a14X4+ ... Fay X,=C1 
aX, as X,+a2X, + aX ... Faz Xx, =C2 


AX Fa,¿X, Fan X5 + aX 4 2... FAm X,=ON 


Eliminando las variables de este sistema debemos construir su matriz 
correspondiente y, añadiendo la columna de los términos independientes, 
nos resultaría: 


a11 212 Aj3 Aj4 -*- Arn cl 
A21 Az Az Az +. Az C2 
Ani An2 An3 Ang +: Ann CN 
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Si se logra triangularizar esta matriz la tendríamos preparada para ob- 
tener la solución del sistema por el método de Gauss-Jordan. Triangulari- 
zar una matriz es un problema que vimos en el tema anterior: aquí el mé- 
todo es semejante, sólo se diferencia en pequeños detalles. Lo primero es 
ver si el sistema es compatible, es decir, si tiene solución. La compatibili- 
dad de un sistema se puede analizar calculando su determinante. Si éste 
es nulo, el sistema es incompatible. Los elementos que componen el de- 
terminante son los coeficientes de las variables a hallar. 

Esquematizando el proceso en el organigrama 9. 


Organigrama 9 


Cada rectángulo representa un proceso diferente; todos ellos encade- 
nados resuelven el problema planteado. 

El programa resultante va a ser especial. Se aplicará el método estruc- 
turado de Bertini para programarlo. Si no está usted habituado a la pro- 
gramación estructurada no se asuste: simplemente este programa va a ser 
aparte de una aplicación matemática, un ejemplo de programación estruc- 
tural. Sin profundizar en el tema, lo esencial del método de Bertini es el 
plantear el problema subdividiéndolo en módulos. De tal forma que cada 
parte o subrutina nos resuelve de manera sencilla una parte del programa 


LL 


principal. Este se compondrá, lógicamente, de varios GOSUB. La lista de 
instrucciones principal sería: 


JA DE COMI 
IMA TE CUOTA 


A PARA DECPETARINOS VALORES DE 3 
A FARA DESFESAR E0OS VALDUREDO DE A 


Los GOSUB 200 y 300 son conocidos por el lector en capítulos anterio- 
res los demás se analizarán uno por uno. Veamos la subrutina de trian- 


Lo primero es convertir el término A(J,J) en 1. Naturalmente, es divi- 
dir por A(J,J), lo cual implica, según las propiedades referidas al princi- 
pio, que esta operación se repetirá a lo largo de toda la fila. 

Una vez conseguida la unidad, se restará a cada fila inferior a J el valor 
de la fila J multiplicada por A(J,J), de tal forma que todos los elementos 
pertenecientes a la columna J, y por debajo de la diagonal, serán nulos. 

Reiterando este proceso a todas las filas la matriz queda triangularizada. 

El proceso de calcular los distintos valores de X es sencillo. Consistirá 
en un bucle que vaya despejando las ecuaciones en proceso inverso, em- 
pezando desde la última hasta terminar con la primera. Las instrucciones 
serían las siguientes: 


La rutina de impresión (GOSUB 600) ya es conocida. Es simplemente 
un FOR que vaya imprimiendo cada uno de los elementos de la matriz X. 

Una vez presentado el método de Gauss-Jordan, pasemos al método de 
Gauss. Al igual que en el anterior la clave era triangularizar la matriz, en 
éste se intenta lograr convertir la matriz de coeficientes en la matriz sin- 
gular. Quedaría: 


X =C1 
X, =C2 
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Por tanto, una vez realizada esta operación resulta X,=C,, X,=C, y así 
sucesivamente. 

Si se comparan los resultados a obtener de uno y otro, se comprueba 
que los planteamientos serán similares. 

Hasta el paso de triangularizar la matriz serán comunes los dos méto- 
dos. A partir de aquí empiezan a diferenciarse. Si en el de Gauss-Jordan se 
despejan las incógnitas, en éste se procederá a triangularizar la matriz 
superiormente. 

De esta forma la matriz quedaría reducida a la diagonal principal, con 
todos sus elementos igual a 1. 

Si se sustituye un proceso por otro, en la operación del ordenador no 
se advertirán diferencias de tiempo. 

El programa del Método de Jordan aparece en el amente listado: 


DDO DE GAUSS-J! 
JO y DADO")! 


Como observará, la explicación dada al método de Gauss-Jordan es apli- 
cable a éste. La diferencia viene dada por la línea 500 y sucesivas. En el 
último método se repite el proceso de triangularización, pero inversamen- 
te; esto se consigue intercambiando los límites de los FOR. 

Con estos dos métodos se pueden resolver los sistemas de N ecuacio- 
nes con N incógnitas. Existe otro llamado de Cramer. Para aplicar este mé- 
todo es necesario: 


e Calcular el determinante (A) de la matriz de coeficientes. 

e Calcular los determinantes (Ai). Estos (Ai) se obtienen reemplazan- 
do en la matriz A la columna i, por el vector de los términos 
independientes. 


Una vez hechos estos cálculos se obtendrían los valores de X,, X,, ..., 
X,, aplicando las siguientes fórmulas: 


x= 141 
|A | En A 
¡As 1 en general = A 
al 


381 


Conseguir que el micro sea capaz de hacer uso de este método no re- 
presenta ninguna dificultad. 

Dejo al lector el conseguirlo: tenga en cuenta que todo se basa en ha- 
llar los distintos determinantes, simplemente habrá que ir jugando con las 
columnas respectivas. 
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COMBINATORIA 42 


ARA el posterior estudio de estadística y probabilidad, son 
necesarias las funciones matemáticas tratadas en este 
capítulo. 

El conocimiento del número de casos posibles de un 
experimento es imprescindible en los cálculos probabilís- 
ticos. Empecemos con las VARIACIONES con REPETI- 
CION. Son necesarios dos datos: 


— Número de elementos: (N). 
— Cantidad de elementos por combinación: (K). 


Sabiendo esto: VRn, k=n*. 


A veces es necesario conocer todas las posibles combinaciones. La ma- 
nera más sencilla es la formación de un árbol. Para ello hacemos: 


aja, 


aja, * a,aza, 


== a,a3a, 


aja, 
7 aja3a, ::* An-1 


a,2,2, 


3833 


¿Cómo controlar este árbol con el ordenador? 

Para programarlo hay que saber qué orden de variación con repetición 
se refiere. El programa siguiente controla el caso de N elementos tomados 
de 4 en 4. 


Si se introdujeran los datos en orden o se utilizara una subrutina de cla- 
sificación al principio, los resultados obtenidos saldrían ordenados. Para 
modificar este programa, de tal forma que sirva para otro tipo de orden 
(K), sólo habría que realizar las siguientes rectificaciones: 


a) Introducir más FOR, hasta que el número de ellos sea igual a K. 
b) En la instrucción 100, añadir tantos A(I) como FOR. 


Puede ocurrir que deseemos eliminar las combinaciones en las que se 
repitan algunos elementos. Estamos hablando de las variaciones ordina- 
rias. En ellas, dos variaciones se diferencian en el orden de colocación de 
los elementos, y en que tienen distintos elementos. 

Basándonos en el programa anterior, sólo habría que añadir las siguien- 
tes instrucciones: 


Con estas nuevas instrucciones se logra que el ordenador salte al si- 
guiente elemento, cuando reconoce dos iguales. Si en las variaciones con 
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repetición URn,k= n*, en las variaciones ordinarias 
Unk=N(n—1)(m—2)...(n—(K —1)). Ver organigrama 10. 


Organigrama 10 


Un caso especial es cuando n=K, estamos entonces en las permutacio- 


nes de n elementos: 
Pn=Un,nn=n(n— 1)(m-2)...3.2.1. 


Al producto 1.2.3...(n—2)(n—1).n se denomina factorial de n(n!). Si se 
quisiera conseguir que el ordenador fuera capaz de calcularlo, simplemen- 
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te se cambiaría la K por N en el programa anterior. Dejo al lector esta mo- 
dificación; no obstante, más adelante aplico esta subrutina en diferentes 
procesos. 

Cuando lo que se desea es elegir un número determinado de elemen- 
tos de un conjunto determinado, estamos hablando de las COMBINACIO- 
NES. Al igual que las variaciones y permutaciones, existen con y sin 
repetición. 

La formación de las distintas combinaciones daría lugar a un árbol del 
siguiente tipo: 

Supongamos A=(a, b, c, d, e) tomado de 3 en 3. 


abc 
LL sl 
ab £ abe 


cede 


386 


Como se deduce, una combinación se distingue de otra en los elemen- 
tos que la integran, no en el orden de colocación. 

Utilizando el programa aplicado a las variaciones, bastaría rectificar el 
valor inicial de los índices de los FOR. Cada FOR empezará con el valor 
de la variable anterior +1, impidiendo de esta forma las posibles repeti- 
ciones. El proceso consecuente se refleja en la siguiente lista de ins- 
trucciones: 


| 
| 


Si se quisiera saber el número de combinaciones, ésta se halla a partir 
de la fórmula: 


E TOS) m)! 


Dejo al lector interesado el estudio del programa para este cálculo. 
Como ayuda indicaré que Cn,k=Cn,n—k, esto es práctico cuando K es 
grande. Se sustituirá K por n—K, ahorrando cálculos innecesarios. 

Si tuviese repeticiones, la fórmula representativa sería: 


CRN,K=Ck+N—1,K 
Utilice el cálculo de Cn,K como subrutina y aplíquelo a CRn,K. 
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Una consecuencia de los números combinatorios es el llamado trián- 
gulo de Tartaglia. Este especial triángulo es de la forma: 


(3) (+) (2) (2) 


Sustituyendo éstos por su valor numérico nos quedaría: 


Lo particular de este triángulo es que cualquier 


ES 
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Es decir, cada elemento es la suma de los dos que están sobre él (a dere- 
cha e izquierda) en la línea superior. Podremos construir un programa ca- 
paz de calcular este triángulo utilizando esta fórmula recursiva. Hay que 
tener en cuenta los siguientes puntos: 


1.2 Determinar hasta qué línea queremos obtener el triángulo. 
n Ñ 
2.2 Tener en cuenta que (n =1 y (al = 1 para cualquier n. 


3.2 Para poder hacer uso de esta fórmula es necesario ir almacenan- 
do los datos anteriores. Para ello guardaremos los resultados en una ma- 
triz NXN, para su uso y posterior impresión del triángulo. 


Fíjese cómo se establecen los límites de los FOR. Debido a que la tabla 
es NXN, habría filas en las que el recorrido total es innecesario; por ello 
se limitan hasta K—1. El resultado para un N cualquiera: 


El binomio de Newton es otra aplicación de los números combinato- 
rios. Este desarrollo utilizado para el cálculo de la potencia M de un bino- 
mio viene dado por: 


(+ay= [9) ar-x0+ []) aro xi+... + [9] -X> 
o bien de manera simplificada: 
(A+Ay3=E (5) Aro xo 
A partir de este término general y aplicando lo aprendido, se logrará 


que el ordenador sepa elevar a una potencia M cualquier binomio. Ten- 
dremos en cuenta las siguientes consideraciones: 


a) Como aplicación concreta utilizaremos un binomio del tipo (X+A). 
b) El cálculo del número combinatorio se hará por medio de una 
subrutina. 


Organigrama 11 
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El programa deducido del organigrama 11: 


Dejo al lector un caso interesante: (A+Bi)? Como sugerencia indicare- 
mos que es similar el proceso, sólo habrá que tener cuidado con el signo. 
(¡Recuerde! que (i)?=1 y (1)?"+!=—1.) 


ESTADISTICA (1) 


UANDO Herman Holletith en 1890 ideó la actualmente su- 
perconocida tarjeta perforada, lo hizo con un fin: realizar 
el censo. Es decir, intentar mecanizar los procesos esta- 
dísticos. Tenía sobradas razones para ello: los cálculos es- 
tadísticos son sencillos pero numerosos. En el momento 
que la recogida de datos es amplia, las operaciones se 
multiplican. 


A lo largo de los dos capítulos siguientes se tratará de 
HR aprovechar la cualidad más importante de nuestro orde- 
nador (la rapidez al operar) para nuestros cálculos estadísticos. Iniciare- 
mos hallando las medidas centrales y presentando las representaciones 
gráficas (capítulo 13) y medidas de dispersión (capítulo 14), concluyendo 
con un análisis estadístico completo. 


En todos los programas se tendrán las siguientes consideraciones: 


e Los datos se almacenarán en una tabla (en realidad se leerían de un 
fichero). De todas formas, los procesos son similares. 


e Todos los cálculos son relativos al espacio muestral. Por tanto, será 
necesario ir acumulando los distintos valores aplicables a estos cálcu- 
los para su posterior uso. 


No se controlan los datos de entrada, se darán por revisados. Esto 
ocurre en la realidad; la verificación de datos suele ser antes de los pro- 
cesos automatizados. 


Comencemos con los cálculos de las frecuencias. Estas no son más que 
el número de veces que se repite determinada variable. 


Para la frecuencia absoluta se van incrementando los distintos acumu- 
ladores según el dato leído. Como no se puede tener indefinidos acumu- 
ladores hace falta limitarlos. Observe el organigrama 12. 
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Organigrama 12 


Estadísticamente se utiliza los intervalos de clase. Si a nuestro ordena- 
dor le introducimos los límites inferiores y superiores y, además, la ampli- 
tud del intervalo, podremos crear una tabla de acumuladores cuya di- 
mensión 


Valor máximo-valor mínimo 


0 Amplitud del intervalo 


El desarrollo para un dato determinado es claro: comparar el valor del 
dato con los extremos inferiores y superiores de cada intervalo; cuando se 
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vas como absolutas. El proceso es idéntico, la única variación es la acu- 
mulación de las frecuencias anteriores a la hallada. Dejo al lector el cálcu- 
lo de estas frecuencias. 


Formación 
Tabla Frecuencias 


Obtención del 
Valor Máximo 


Organigrama 13 
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Comencemos con las medidas llamadas centrales. Estas son: 


e Media aritmética. 
e Moda. 
e Mediana. 


La moda m, de una serie estadística, es el valor de mayor frecuencia. 

Para su estudio es necesario dos subrutinas ya vistas: una encargada del 
cálculo de las frecuencias, vista en este capítulo, y otra de clasificación 
para obtener la frecuencia máxima (capítulo 1). Estructuraremos nuestro 
programa de acuerdo al organigrama 13. 


El listado resultante y la solución para un caso concreto: 


mn mm” LA MM 
LO DE LA MOI] 


La media aritmética viene dada por la fórmula: 


k 
MA = + 2 Xi-fi 
en donde Xi son los distintos valores de X, fi su frecuencia y N el número 
total de valores, lógicamente N=fi+f2+ ... +fk 
Nos encontramos de nuevo con el esquema de trabajo indicado al prin- 
cipio del capítulo. El programa es similar al de la moda, la variación será 
la sustitución de la subrutina de clasificación por una que calcule el algo- 
ritmo de la media. Aprovechando la lista de instrucciones anterior e intro- 
duciendo nuestra modificación, el resultado sería: 


11 


Queda por ver la Mediana. Su definición estadística es: 
Se llama mediana de una serie estadística al valor M, tal que existe igual 
número de observaciones mayores que menores que M. 
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Si se observa, no hay más que coger el elemento N/2 de nuestra matriz 
de datos. Pero si no queremos un valor erróneo, es necesaria una clasifi- 
cación previa de esta tabla. 


El proceso reflejado en el organigrama 14 será más fácil de asimilar: 


La resultante 
de esta clasificación 
es la tabla C(N) 


Organigrama 14 


Conocemos, por el capítulo 1, cómo efectuar una clasificación (SORT) 
de un conjunto determinado. Si se aplica esta subrutina a nuestros datos 
iniciales, la obtención de la mediana no tiene complicación alguna. Antes 
de escribir nuestra lista de instrucciones, se debe tomar en consideración 
un detalle importante. Pese a la rapidez de operación del ordenador, es ne- 
cesario el ahorro de tiempo, máximo cuando los datos pueden ser nume- 
rosos. Se necesita el elemento N/2 de la tabla ya ordenada; si se clasifica 
hasta dicho elemento ahorraremos la clasificación de los términos poste- 
riores. Así en lugar de N pasadas, bastará N/2 con la consiguiente ganan- 
cia de tiempo. 


70 NEXT J 

80 REM CLASIFICAFION DE LOS DATOS 

90 FOR J=1 TO N 

100 FOR K=J TO N 

110 IF D(J)<D(K) THEN LET X=D(J):LET-D(J)=D(K):LET 
D(K)=X 

120 NEXT K 

130 NEXT J 

140 LET S=INT (N/2)+1 

150 PRINT "LA MEDIANA ES: ";D(S) 

999 END 


A lo largo del capítulo se han visto conceptos estadísticos de las llama- 
das medidas centrales. A continuación trataremos las representaciones 
gráficas más usuales. 

Existen varios tipos de gráficas estadísticas; trataremos dos de ellas: la 
de barras y la de sectores. 

La de barras consiste en representar los dos ejes y sobre el eje X elevar 
barras cuya elevación sea el porcentaje de la variable. El programa resul- 
tante será: 


10 REM GRAFICAS DE BARRAS PARA EL AMSTRAD 

20 MOVE 0,0 

30 INPUT "NUM. DE ELEMENTOS :"¿N 

40 IF N>200 OR N<i THEN PRINT "ERROR: ENTRE i Y 20 
0":G0TO 30 

SO DIM TABLA (N) 

60 INPUT "TECLEE DIFERENCIA ENTRE ESCALAS: "¿V 
70 FOR I=i TO N 

80 PRINT “DATO "1: INPUT TABLA(I) 

90 IF TABLA(I)>25%*V THEN PRINT "ERROR: VALOR SUFPERI 
OR A"325*V:GOTO 80 

100 NEXT 1 

110 CLS 

120 BARRADO=INT (200/N) 

140 PLOT 100,200: DRAW 100,400 

150 FOR I=BARRADO TO 200 STEP BARRADO 

160 PLOT 100+1,195:DRAW 100+1,200 

LOA 

180 FOR J=V TO 200 STEP BARRADO 

190 PLOT 93,200+J:DRAW 100+1,200 

200 NEXT J 

210 J=0 

220 FOR I=BARRADO TO 200 STEP BARRADO 


100 


230 J=J+1 

240 CIMA=25*V 

250 PORC=INT(TABLA (J)*100/CIMA) 
260 SUBIDA=INT (PORC*CIMA/100) 
270 PLOT 100+1,200 

280 DRAW 100+1,200+SUBIDA 

290 NEXT I 


Otro tipo de gráficas consiste en dividir la circunferencia en secciones, 
cuyos ángulos son proporcionales de las variables. El razonamiento es si- 
milar al anterior. 


10 REM GRAFICA DE SECTORES PARA AMSTRAD 
20 INPUT “CUANTOS ELEMENTOS: "¿N 

30 DIM TABLA(N) 

40 FOR I=i TO N 

50 PRINT "ELEMENTO “51: INPUT TABLA(I) 
55 NUM=NUM+TABLA (1) 

60 NEXT I 

7O CLS 

80 DEG 

90 ORIGIN 320,200 

100 FOR I=1 TO 360 

110 PLOT 100*C0S(I),100*SIN(I) 

120 NEXT 1 

130 DIM GRAD(N) 

140 FOR I=1 TO N 

150 PORC=INT(TABLA(1)*100/NUM) 

160 GRAD(1)=INT (360*PORC/100) 

170 NEXT 1 

200 G=0 

210 FOR J=1 TO N - 

215 ORIGIN 320,200 

220 DRAW 100*COS(GRAD(J)+6) , 100*SIN(GRAD (J) +6) 
230 G=G+GRAD (J) 

240 NEXT J 


Existen otros tipos de gráficas; su mayor o menor dificultad depende- 
rán de la pericia del lector. ¡Animo! 
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ESTADISTICA (ID) 


A 


MEDIDAS DE DISPERSION 


N el capítulo anterior se trató de las medidas centrales y, 
posteriormente, se logró hacer gráficas estadísticas con el 
ordenador. Ahora se verán las llamadas medidas de dis- 
persión. 

La utilidad de éstas viene dada por la necesidad de 
completar la información conseguida por las medidas 
centrales. Si se observa la siguiente tabla: 


Calificaciones Calificaciones 


TETERA 


Media X = 2,90 Y = 2,85 


Se ve cómo las medias son casi iguales, pero en la tabla B, los datos es- 
tán más dispersos. Por esta razón se pensó en la utilización de otras medi- 
das, que informarán sobre esta dispersión. Las medidas de dispersión son 
siempre respecto a la media aritmética. Ya tenemos un punto de apoyo a 
la hora de estructurar los programas. Habrá que echar mano de los pro- 
cesos y subrutinas del tema anterior para el cálculo de las medidas de 
dispersión. 
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La primera en tratar es la desviación media. Para ello, vamos a definir 
el concepto de Recorrido: 


— Recorrido de una serie estadística es la diferencia entre el valor 
máximo y el valor mínimo. 


Fíjese en el organigrama 15. 


Organigrama 15 


Tanto la recogida de datos como su posterior clasificación no deben 
ofrecer dificultad alguna, ya que ha sido tratado en temas anteriores. Es- 
pero que el lector no encuentre ningún obstáculo en la confección del 
correspondiente programa. Lógicamente UMax=A(1) y Umin=A(N), siem- 
pre y cuando esté ordenada la tabla de mayor a menor; en caso contrario, 
la asignación del valor máximo y mínimo sería la inversa. 

Una vez conocido el recorrido continuemos con el concepto de desvia- 
ción media. 

Al estudiar la media aritmética, se vio que era una medida de tenden- 
cia central; pues bien, la desviación media indica el valor en que los dife- 
rentes elementos de una serie se alejan de este punto central. 

Lógicamente lo ideal es que esta medida fuese igual a 0, pero esto no 
sucede nunca. 

Cuando el espacio muestral es pequeño, la fórmula estadística más 
apropiada es: 


y | Xx-Xm| 


DM = N 
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en la que: 


X  = Valor de los distintos elementos. 
XM = Valor de la media aritmética. 
N  = Número de elementos. 


Reflejemos el desarrollo en el organigrama 16. 


Organigrama 16. 


El programa deducido del organigrama 16. 


10 REM CALCULO DE LA DESVIACION MEDIA 
- 20 REM ENTRADA DE DATOS 


TOS 


Cuando el número de datos es numeroso sería conveniente dividir el 
espacio muestral en intervalos. Esto ya lo hicimos en el tema anterior; aho- 
ra bien, el cálculo de la desviación media varía. El camino a seguir sería: 


1.2 Se distribuye la serie por clases. 

2.2 Se hallan las frecuencias (X), el número de elementos que perte- 
necen a cada clase. 

3.2 Se halla el punto medio, marca de clase, de cada intervalo. 

4.2 Se halla la media. 

5.2 Se calcula (X—XM), X=Marca de clase. 

6.2 Se multiplican las frecuencias por (X—Xm). 


La fórmula nos queda: 


2f- | xX—-Xm| 
N 


Fíjese en el organigrama 17 e intente construir el programa. 

Todos los procesos son conocidos, luego no creo que presente ningu- 
na dificultad la realización de la lista de instrucciones. 

Otra medida de dispersión para hallar la homogeneidad más o menos 
acentuada de un conjunto de elementos es la desviación típica o estándar. 

La fórmula es la siguiente: 
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Organigrama 17 


en la que 


Y | X-Xm |? = Suma de los cuadrados de las desviaciones. 
N = Número de elementos. 


Ya que en la desviación media para el caso de pocos elementos, en ésta 
lo haremos para un gran espacio muestral. 
En este caso la fórmula quedaría: 


donde: 


X = Punto medio de cada intervalo. 
F = Frecuencia de cada clase. 
N y Xm son el número de elementos y la media aritmética. 
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Los pasos son los siguientes: 


1.2 Recogida de datos. 

2.2 Hallar las frecuencias de cada clase. 

3.2 Calcular y almacenar en una tabla las marcas de cada clase. 
4.2 Aplicar la subrutina para hallar la media aritmética. 

5.2 Hallar f- | X—Xm]? e ir acumulando estos colores. 

6.2 Aplicar la fórmula 6 e imprimir el resultado. 


Todos los pasos han sido calculados anteriormente, sólo el punto 5." es 
diferente, pero es similar al estudiado en el programa anterior. 

Fíjese en el programa; los comentarios (REM) irán indicando cada uno 
de los diferentes pasos a realizar: 


10 REM CALCULO DE LA DESVIACION TIFICA 

20 INPUT "NUMERO DE DATOS";¿N 

25 LET SU=0 

30 DIM A(N)> 

40 PRINT "INTRODUZCA DATOS: " 

50 FOR J=i TO N 

60 INPUT D(J): SU=SU+D(J) 

70 NEXT J 

735 LET ME=SU/N 

80 INPUT "INTRODUZCA AMPLITUD DE LOS INTERVALOS: "; 


90 INPUT "INTRODUZCA LOS VALORES MAXIMO Y MINIMO: " 
¡ MX, MN 

100 LET L=INT(MX-MN)/A 

110 DIM C(L) 

120 FOR J=1 TO MN 

130 FOR K=1 TO L 

140 IF D(J)<(MN+Kx*A) THEN LET C(K)=C(K)+1:GOTO 160 
150 NEXT K 

160 NEXT J 

170 REM FORMACION DE LA TABLA DE LOS PUNTOS MEDIOS 
180 DIM M(L) 

190 FOR J=1 TO L 

200 LET M(J)=MN+A* (J-1)+A/2 

210 NEXT J 

220 REM CALCULO DEL CUADRADO DE LAS DESVIACIONES 
230 FOR J=1 TO L 

240 LET D=ABS(M(J)-ME) 

250 LET SU=SU+D"2*C (J) 

260 NEXT J 

270 LET DE=SQR(D/N) 

280 PRINT"LA DESVIACION TIPICA ES: "DE 


1083 


Un método estadístico para ahorrar operaciones es aplicar el método 
abreviado de elegir una media supuesta. A partir de esta media calcular 
las desviaciones (D) y aplicar la fórmula siguiente: 


en la que 


I = Intervalo de clase. 
N = Número de elementos. 
Simplemente es modificar el anterior programa de la siguiente manera: 


1.2 Pedir la supuesta media; por tanto, eliminar la subrutina de la 
media. 
2.2 Sustituir esta fórmula donde se aplique la otra. 


Dejo al lector comprobar la rentabilidad de un método y de otro. Hay 
que tener en cuenta que siempre es necesario probar los métodos para va- 
rios grupos de elementos. Muchas veces el empleo de distintos métodos 
en procesos informáticos varían según la cantidad de datos a tratar. Un 
ejemplo claro son los distintos métodos de clasificación; su elección de- 
pende de variadas fórmulas que midan la rentabilidad del empleo de uno 
u otro, en relación al número y tipos de datos. 

A veces es interesante conocer la desviación en porcentaje. Habría que 
añadir a nuestro programa los siguientes puntos: 


— Calcular Ma + o. 

— Ver cuántos elementos están en (Ma—0, Ma+0). 

— Calcular el porcentaje que representan estos elementos en el espa- 
cio muestral. 


Traducido a instrucciones sería: 


300 REM SUBRUTINA PARA EL CACULO DE PORCENTAJES 
310 FOR J=1 TO N 

320 IF D(J)>(ME-DE) AND D(J)< (ME+DE) THEN LET CN=C 
N+1 

330 NEXT J 

340 PRINT "EL PORCENTAJE ES: ";CN*100/N 

350 RETURN ] 


Después de lo visto en estos dos capítulos tenemos todos los elementos 
necesarios para un análisis completo estadístico. 
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Supongamos que tenemos los siguientes datos: 


APLICANDO LOS DISTINTOS PROGRAMAS: 
LA MODA ES:170.8 

LA MEDIA ES:170.7 

LA MEDIANA ES:170.7 

LA DESVIACION TIPICA ES:10 

EL FORCENTAJE ES 297 


ESPACIO A” > 


L estudio de los vectores y su posterior aplicación en la Fí- 
sica hace necesario que se trate en este libro. 

Un vector es un segmento orientado, de tal forma que 
es necesario conocer la dirección, módulo y sentido del 
mismo. 

Lo primero es saber operar con vectores; en el siguien- 
te programa se aprenderá a realizar estos procesos: 


1.2 Cálculo del módulo. 

2.2 Suma de vectores. 

3.2 Producto de un vector por un escalar. 
4.2 Producto escalar. 


Para este último punto se calculará, previamente, el valor del ángulo 
formado por dos vectores. La subrutina correspondiente es: 


| 


ONDIVNNAD” 
JIOIKCAUAZ 2 (1123 SORK 
LINA 


SO PRINT "3-PRODUCTO DE UN VECTOR POR UN ESCALAR" 
60 PRINT "4-PRODUCTO ESCALAR" - 

70 INPUT OP:IF OP<1 OR OP>4 THEN GOTO 70 

80 IF OP=1 THEN GOTO 100 

82 IF OP=2 THEN GOTO 200 

84 IF OP=3 THEN GOTO 300 

86 IF OP=4 THEN GOTO 400 

90 END 

100 INPUT "INTRODUZCA COORDENADAS DE LOS EXTREMOS: 
*7XO,YO,X1,Y1 

110 LET MO=SQR (X1-X0)”2+SQR (Y1-Y0)”2) = 
120 PRINT "EL MODULO ES";MO == 
130 RETURN 

200 INPUT "INTRODUZCA LAS COORDENADAS DE LOS VECTO 

RES: "¿V1,U1,V2,U2 === 

210 PRINT "EL VECTOR SUMA ES!";V1+V23"y"¿U1i+U2-————— 
220 RETURN = : = - 

300 INPUT "INTRODUZCA LAS COORDENADAS DEL VECTOR Y 
EL ESCALAR"yV1,U1,A z 

310 PRINT "EL RESULTADO ES: "¡V1x*A;U1x*A ES 

320 RETURN : AAZRAA 
400 INPUT "INTRODUCIR LAS COORDENADAS DE LOS DOS Y 

ECTORES: ";V1,U1,42,U2 

410 LET M1=SQR(V1"2+U172) 

420 LET M2=SQR(V2"2+U2"2) 

430 REM LLAMADA A LA SUBRUTINA DEL CALCULO DEL COS 

ENO INDICADA ANTES 

440 GOSUB 500 

450 LET PR=M1x*M2x*CS 

460 PRINT "EL PRODUCTO VECTORIAL ES: ";PR 

470 RETURN 


Como se ve, el programa principal (10-90) se compone de un menú 
para la elección de operación. Según ésta serán necesarios unos datos ini- 
ciales, introducidos por medio de los input's. Las operaciones en sí no pre- 
sentan ninguna dificultad, sólo hay que reseñar la llamada a la subrutina 
del coseno, vista al principio del programa, necesaria para el cálculo del 
producto escalar. 

Una vez que hemos aprendido a operar con vectores, dejaremos para 
el siguiente capítulo los problemas relacionados con bases y cambio de las 
mismas. Ahora se trata de aprovechar al máximo las coordenadas de un 
punto. 

Dados N puntos, con sus respectivas coordenadas es posible hallar el 
área del recinto limitado por ellos. Esto es gracias a la fórmula de Gauss. 
Esta es: 


AMEX E IN) AO Y dED) 
Z 


AREA = 
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La única limitación es la imposibilidad de que se corten los lados, de 
tal forma que el perímetro se cruzase. 
Ver organigrama 18. 


Organigrama 18 


Una vez conocida la fórmula para N lados nos centraremos en una fi- 
gura ya estudiada en este libro: el triángulo. 

Conocidos tres puntos, se trata de calcular todos los puntos notables 
del triángulo, tales como: baricentro, incentro, ... 

El procedimiento consistirá en ir aplicando las distintas propiedades 
ya vistas. Todas las fórmulas empleadas son aplicación directa de fórmu- 
las matemáticas. 
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A B 


El programa nos resultaría: 


Los datos iniciales (20-40) serán, lógicamente, las coordenadas de los 
tres puntos. Los resultados se obtienen aplicando las fórmulas geométri- 
cas correspondientes. Partiendo del cálculo de la longitud de los lados (dis- 
tancia entre dos puntos) se obtienen el baricentro y el incentro (100-130). 
Gracias al cálculo de los ángulos es posible saber el circuncentro y orto- 
centro (200-260). 
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LA RECTA 6 


N el capítulo 8 tratamos las ecuaciones de primer grado, 
posteriormente (en el capítulo anterior) vimos las coor- 
denadas en el espacio vectorial. En éste, uniremos ambos 
conceptos para tratar la línea recta. 

Intentaremos lograr que nuestro ordenador maneje 
conceptos como paralelismo, perpendicularidad, etc. 
Como habitualmente hacemos, comenzaremos con lo más 
sencillo; tengamos en cuenta que nuestro micro no parte 
con los conocimientos básicos necesarios. 


Hay cinco formas comunes de representar a la línea: 


1.2 Ecuación general: 
AX+by=C 


2.2 Forma explícita: 
Y—Y1=M(X-X,) 


3.2 Paramétricas: 


X=X, +T*V, 
Y=Y, +T*V, 
4.2 Continua: 
XX; mn YY; 
V, V, 


5.2 Función segmentos que determina sobre los ejes: 


> A: A 
E. a" 
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Lo 


primero es conseguir saber traspasar la representación de la recta 


de una a otra. Para ello, según los datos conocidos, pueden darse los si- 
guientes casos: 


qe 
7 
Se 
4, 
2 


Conocidos A, B, C. 

Conocidos dos puntos (X1, Y1) y (X2, Y2) por los que pasa. 
Conocidos un punto (X1, Y1) y el vector dirección (V1, V2). 
Conocidos un punto (X1, Y1) y la pendiente M. 

Conocidos los segmentos sobre los ejes. 


Estas son las posibles opciones para determinar una recta. 

Confeccionaremos un menú para que, según la opción tecleada, nues- 
tro ordenador reclame los datos necesarios. 

En todos los casos hay que conseguir los mismos parámetros A, B, C, 
V1, V2, M, F, G, X1 y X2; según tengamos unos u otros aplicaremos las re- 
laciones matemáticas oportunas. He aquí el programa; las anotaciones de 
la derecha te servirán para seguir el desarrollo matemático: 


10 REM PROGRAMA PARA EL TRASPASO DE UNA FORMULA DE LA RECTA A OTRA 


-20 PRINT "ELIJA OPCION: " 
30 PRINT "1-CONOCIDO A.B,C" 
40 PRINT "2-CONOCIDO DOS PUNTOS" 
50 PRINT "3-CONOCIDO UN PUNTO Y EL VECTOR DIRECCION" 
- 60 PRINT "4-CONOCIDO UN PUNTO Y LA PENDIENTE" 
70 PRINT "5-CONOCIDO LOS SEGMENTOS SOBRE LOS EJES" 
80 INPUT OP: IF OP<1 OR OP>6 THEN GOTO 80 
- 90 1F OP=1 THEN GOSUB 200 
91 IF OP=2 THEN GOSUB 300 
92 1F OP=3 THEN GOSUB 400 
93 IF OP=4 THEN GOSUB 500 
95 1F OP=5 THEN GOSUB 600 


GOSUB 700 

END 

INPUT "A,E,C: "3A,B,C 

LET M=-B/A 
FLET X1=1:LET Yi=(AXX1+C)/-B 


LET X2=2:LET Y2=(AkX2+€)/-B 
LET Vi=X2-X1:LET V2=Y2-Y1 
LET F=C/A:LET G6=C/B 


RETURN 


INPUT "INTRODUZCA X1,Y1,X2,Y1 "3X1,Y1,X2,Y2 
LET Vi=X2-X1:LET V2=Y2-Y1 

LET M=-(Y2-V1)/(X2-X1) 

LET A=Y2k(XZ8X1)-X1X(Y2-Y1) 


RETURN 


INPUT "INTRODUZCA X1,Y1,V1,V2: "3X1,Y1,V1,V2 
LET M=-42/V1 

LET A=V2:+LET B=V1 

LET C=Y1XV1-X1XV2 

LET F=C/A:LET G=C/B 

RETURN 

INPUT "INTRODUZCA X1,Y1,M "3X1,Y1,M 

LET X2=2:LET Y2=2 


Ty 


515 LET V2=-V1XM 
520 LET A=Y2:LET B=-V1 
530 LET C=Y1RV1-X1XV2 
540 LET F=C/A:LET G=C/B 
550 RETURN 
600 INFUT "INTRODUZCA LOS DOS SEGMENTOS F,G ";F,G 
610 LET C=F:LET A=1:LET B=C/G 
620 LET M=-B/A 
630 LET Xi=1:LET Y1=(AKX1+C)/-B 
640 LET X2=2:LET Y2=(AXX2+C)/-B 
650 LET Vi=X2-X1:LET V2=Y2-Y1 
660 RETURN 
700 REM IMPRESION 
710 PRINT "EN FORMA GENERAL ES: " 
720 PRINT Aj "X+"Bj"Y+"303"=0" 
730 PRINT "EN FORMA EXPLICITA ES: " 
790 PRINT "Y="7Y13"="3M3"(X="3X13 00" 
750 PRINT "EN FORMA CONTINUA ES: " 
760 PRINT "X-"3X13"/"3V13"= V-"3Y13"/"3V2 
770 PRINT "EN FORMA PARAMETRICA ES: " = 
780 PRINT "X="3X13"+ Tk"3V1 
790 PRINT "Y="5Y13"+ TX"3V2 
800 PRINT "X/"3F3"+Y/"3G3"=1" 
810 PRINT "X/"3FpM+Y/"3G3"=1" 
820 RETURN 


Por medio del menú principal se bifurca a la opción deseada. En ésta 
se pedirá los datos necesarios y a partir de éstos se calculan el resto de los 
datos. En todas las direcciones dadas por el GOSUB se sigue este proceso; 
las fórmulas se obtienen despejando el dato, hallar de una u otra forma de 
representación, siempre teniendo en cuenta los datos iniciales. 

Obtenida la recta, juguemos con ella; o, mejor dicho, consigamos que 
nuestro ordenador lo haga. Comencemos con la relación de paralelismo y 
perpendicularidad de dos rectas. Para conocer estas relaciones es esencial 
saber el ángulo que forman. El ángulo formado por una recta y el eje OX 
viene dado por la pendiente M(M=tanga). Si deseo una recta paralela, ésta 
deberá tener la misma pendiente. No habría ningún problema en hallarla. 
Con la perpendicular habría algunas dificultades, ya que habría que au- 
mentar (o disminuir) en M/2 radiones. 

Matemáticamente sabemos la condición de que una línea sea perpen- 
dicular a otra, es: 


Pongámoselo difícil a nuestro ordenador y programémosle para con- 
seguir la pendiente progresivamente sin aplicar directamente la fórmula. 
Supongamos conocidos A, B y C inicialmente (si éstos no fueron los da- 
tos, aplicando el programa anterior, los calcularíamos) y nos pidieran 
calcular la paralela y la perpendicular que pase por el punto (Xo, Yo). Con 


az 


la paralela no habría problema y quedará solucionado entre las líneas 70 
y 120. Pero la perpendicular la vamos a hallar de diferente manera. Ire- 
mos girando la recta paralela obtenida, de 10 en 10 grados, hasta conse- 
guir 90%. En este momento, comprobaremos dos cosas: 


1.2 Que la recta obtenida tiene su pendiente igual 
1 
M2 
2.2 La distancia entre el punto (Yo, Yo) y la recta es mínima. 
Para este segundo punto seguiremos los siguientes pasos: 


a) Calculemos el punto de corte (X1, Y1) entre la recta dato y las di- 
ferentes rectas que voy hallando. 
b) Calculo la distancia (X1, Y1) y (Xo, Yo). 


Como comprobación adicional si en vez de parar el giro en 90% se con- 
tinúa hasta 100%, observará cómo empieza a aumentar la distancia. 


MDI 
DO 


—. ITZO LET BD=ecoarRít (Y o 
2 Y LEI DESTA LALA? CTA a 


El cálculo de la paralela (líneas 50, 60) no representa ninguna dificul- 
tad al aplicar directamente las fórmulas de paralelismo. Con el FOR (línea 
90-120) se produce el proceso de girar la línea 100 grados, de 10 en 10. 
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calculando la distancia (98) a la vez que la pendiente (95). Para saber la 
solución se utiliza la instrucción 110. 

Si se atreve, la versión gráfica de este programa puede serle interesan- 
te. Sólo tiene que ir tomando como referencia el punto de corte y lanzar 
líneas (por ejemplo, DRAW en el Amstrad) desde (Yo, Yo) a los distintos 
puntos de corte. 


CURVAS Pe y 


L estudio de las cónicas con sus relaciones es el tema a tra- 
tar en este capítulo. Empezaremos con la más conocida 
de todas ellas: la circunferencia. 

El lugar geométrico de los puntos del plano que equi- 
distan de un punto fijo llamado centro es la circun- 
ferencia. 

¿Qué es un lugar geométrico? 
=== Es el conjunto de puntos que poseen una determinada 
===> propiedad (y además son los únicos puntos que la poseen). 

Un dato a tener en cuenta es el radio. Definido como la distancia del 
centro a cualquier punto de la circunferencia. 


(%o» Yo) 


El radio vendrá dado por la expresión: r=y2+Y?. Elevando al cuadra- 
do se obtendrá la fórmula general de la circunferencia. 


R2=X? +Y? 
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El siguiente programa calculará los valores de los puntos de la circun- 
ferencia. Como todos no se pueden hallar (son infinitos), impondremos un 


intervalo 5/100. 


- O REM CALCULO DE LA TABLA. DE VALORES. DE UNA CIRCU 
— NFERENCIA ===> 

20 INPUT “INTRODUZCA RADIO: AR 
== 
APRISTA EZ === 
59 LEY DE=2*3. 14159/100 >= == 


Ya hemos conseguido hallar una tabla de valores. El dibujar esta cir- 
cunferencia lo dejo al lector. Tenga en cuenta que en los Manuales de la 
mayoría de los micros vienen unas subrutinas de ejemplos; una de ellas 
suele ser el dibujo de la circunferencia. 


La circunferencia puede venir dada por la fórmula llamada general: 


AX?+BY?+ CXY+DX+EY+F=0 


Las condiciones que deben cumplir estos coeficientes son: 


A=1.K D=-—2ak 
B=1.K E=—2bk 
C=0.K F=(a2+b2—r?) K 


Si despejamos a, b y r nos queda: 


Do ja E ju WIPE4AR 


2A 2A 2A 


Es decir, será una circunferencia de centro 
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Dada una recta puede ser interesante saber su posición respecto a una 
circunferencia. Fíjese en la figura y donde pueden verse los tres posibles 
casos: 


Ver organigrama 19. 


Organigrama 19 


123 


La lista de instrucciones no revistirá ninguna dificultad. Simplemente 
habrá que resolver el sistema dado por las ecuaciones de las dos figuras 
geométricas: recta y circunferencia. Nos quedará una ecuación de segun- 
do grado, y analizando su discriminante se obtendrá la posición de la recta. 


10 REM CALCULO DE LA POSICION DE UNA RECTA EN UNA 
CIRCUNFERENCIA 
20 INPUT "INTRODUZCA EL RADIO 3: "¿R 

30 INFUT "INTRODUZCA COEF ICIENTES. DE LA RECTA (A,B yn 
10) 1"594,+B,C - 
40 LET DE=2*3.14159/100 = 
SO REM X“2+ (=-C-AX)“2=R"2. DESARROLLANDO QUEDA: (A” 
2+1)7X+2C0X+ (0*2-R"2) 

60 LET B=2x*C:LET A=A*2+1:LET C=C"2-R”"2 
70 LET DI=B"2-4x*Ax*C : 

= 90 IF DI¿0O THEN PRINT "LA RECTA ES EXTERIOR." 
100 IF DI=0 THEN PRINT"LA RECTA ES TANGENTE" 
110 IF DI>0 THEN PRINT "LA RECTA ES SECANTE" 

999 END 


13X 


Un concepto importante en el estudio de las circunferencias es la po- 
tencia de un punto respecto de una circunferencia. Dado un punto 
P(Xo, Yo), si de él se trazan secantes a la circunferencia se cumple: 
PA.PB=constante, siendo PA el segmento que une P y el primer punto de 
corte, y PB lo mismo con el segundo punto de intersección de la 
circunferencia. 

A este producto se llama potencia de un respecto a una circunferencia. 
Si logramos trazar un secante tal como la figura: 


Se ve que P=d*—r?, ya que PA=d—r y PB=d+r siendo d la distancia 
del punto a la recta, sabiendo que d?=(Xo—a)?+(Yo—b)?, la potencia será: 


P¿,=(Xo—a)?+(Yo—b)?— r? 
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Dadas dos circunferencias se llama eje radical a la recta compuesta por 
los puntos que tienen la misma potencia respecto a las dos circunferen- 
cias. Dado los centros y radios de dos circunferencias calcularemos su eje 
radical al siguiente programa: 


10 REM CALCULO DEL EJE RADICAL DE DOS CIRCUNFERENC 


IAS 
20 INPUT "INTRODUZCA CIRCUNFERENCIA 1 EN FORMA GEN 


ERAL 3 ";D1,E1,F1 
30 INPUT "INTRODUZCA CIRCUNFERENCIA 2 EN FORMA GEN 


ERAL 3: "5D2,E2,F2 

40 PRINT "EL EJE RADICAL DE LAS DOS CIRCUNFERENCIA 
S3-E91" 

50 PRINT “COEFICIENTE DE LA X: "¿D1-D2 

60 PRINT "COEFICIENTE DE LA Y: "¡El-E2- 

70 PRINT "TERMINO INDEPENDIENTE: "¡Fi-F2 

999 END 


Otra cónica conocida es la elipse. Se define como el lugar geométrico 
de los puntos, tales que la suma de las distancias a dos puntos fijos llama- 
dos focos en constante. 


Elementos de la elipse: 


F y F: Focos de la elipse. 

Eje focal: Eje que pasa por los focos. 

Eje secundario: Eje mediatriz del segmento FF. 

Vértices: Puntos donde cortan los ejes a la elipse. 

Radios vectores: Son los segmentos que unen los focos con un punto 


cualquiera de la elipse. 


12.5 


La fórmula de la elipse vendrá dada por: 


e O Y e O 
a? b? a 


donde Xo y Yo son las coordenadas del centro. 


Veamos cómo podemos imprimir una elipse. Para ello reduzcamos la 
fórmula a: 


Hay que tener en cuenta que en la impresora, al menos las usuales, no pue- 
den subir y bajar. El gráfico habrá que hacerlo de tal forma que vayan sa- 
liendo las líneas horizontalmente. Es decir, la elipse quedará como un ba- 
lón de rugby. Al aplicar la fórmula, ésta nos da dos posibles valores de y: 


= AA 
40 PRINT " - 
50 LET. 1=B/20. : - 
60 FOR J=1 TO 21 
65 LET X=IR(J-1)- : 
FO LET V1=(C=(BO2ZAAACZII SAC 
40 ET 12 (CABTZIOS ZOAACZ 
A == : 
566 NEO E == 7 - 


AAA = 


Partiendo desde X=A hasta X=0 tendremos media elipse dibujada. Si 
lo hiciéramos desde X=— a hasta X=a obtendríamos la elipse entera. Se 
va a hacer sólo media elipse, para después repetir el FOR de 0 a A, consi- 
guiendo la otra media. 

El programa no tiene dificultad, bastará controlar los blancos antes y 
después de imprimir los asteriscos. 

¿Dónde estarán los ejes? 

Debido al problema de la impresora, los ejes quedarán invertidos. 

Otro punto a tener en cuenta es la limitación de nuestra gráfica. Pode- 
mos disponer de un ancho de 80 columnas, luego el valor máximo de 
Y será 39 (80— 2/2). Esta cota habrá que tenerla en cuenta a la hora de sa- 
ber el valor de cada posición de la impresora. 

¿Cómo calcularla? 


LZ6 


Muy sencillo, el valor máximo de Y es cuando X=0, sustituyendo este 
valor en la fórmula: Y = 

Si como mucho el ancho puede ser B, cada posición representará un 
valor de: (39/B). Tampoco sirven los decimales, luego habrá que utilizar 
la función INT. Esto originará pequeños relieves en nuestra gráfica. El pro- 
grama quedaría: 


Ya sabe el ordenador manejar circunferencias y elipses, completos 
su formación con la hipérbola. Esta viene representada por la figura: 


1.27 


Su definición es: lugar geométrico de todos los puntos del plano tales 
que la diferencia de las distancias a dos puntos fijos llamados focos es una 
cantidad constante. 


Elementos significativos: 
Focos K y F', son los puntos fijos de la definición. 


Distancia focal: FF. 
Radios vectores: PF y PF”. 
Eje focal: recta que pasa por los focos. 
Se cumple que para cualquier punto P: 
Su fórmula general: 
(X—Xo)? (Y—-Yoy? 
A 
Siendo (Xo, Yo) los coordenados del centro de la hipérbola. 


Al igual que se hizo con las otras dos cónicas se programará al ordena- 
dor para obtener la fórmula según unos datos dados. 


==] 


10 REM CALCULO DE LA FORMULA GENERAL DE LA HIPERBO 
20 INPUT "INTRODUZCA COORDENADAS DEL CENTRO (X, Ya 


SANS “INTRODUZCA VALORES DEA Y Bi 130,5 AS - 
40 PRINT "EL TERMINO DE X"2 ES: “yBoz == 
50 PRINT "EL TERMINO DE Y*2 ES: Az == 
60 PRINT "EL TERMINO DE X ES 2 "3-(XKB"2)- 
70 PRINT "EL TERMINO DE Y ES 2 "¿YRA“2 z 
80 PRINT "EL TERMINO INDEPENDIENTE ES : "3 (E92xX "2 
)(AC2ZRBOZ)(AMZAV OZ) 


Habiendo comprendido la circunferencia y la elipse no entraña dificul- 
tad alguna el programa anterior. 

Un caso interesante es cuando dan una ecuación y piden decir qué cla- 
se de cónica es. 

Para ello, se basará el razonamiento en la siguiente propiedad ma- 
temática: 


Sea AX?+BXY+CY* FDX+EY+F 
Se hace: A*C*F— 2*B*E—A*E?—C*D?—F*B*=DETER 


Si DETER >0 es hipérbola 
DETER <O0 es elipse 
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El programa deducido será: 


DULLA DATO: 


YTIIDOITIEE LAT ECMINO EIA AEE A cd 
O A er 


O INPUT "EL TER 


Queda para el lector el intentar imprimir la hipérbola y la circunferen- 
cia, el proceso es similar al de la elipse. Después de este capítulo nuestro 
ordenador domina la geometría... o al menos es menos torpe de lo habitual. 
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ESTUDIO DE poda 


L estudio de las gráficas de las funciones es cuestión de ex- 
periencia y suerte. La primera es imprescindible para apli- 
car a simple vista los conocimientos teóricos en los pun- 
tos adecuados, evitando, de esta manera, pruebas con va- 
lores innecesarios. 

Suerte, sobre todo al estudiar la continuidad. En deter- 
minadas funciones la discontinuidad está clara; este es el 
caso, por ejemplo, de las funciones con denominadores. 
Hallando las raíces de éstos, ya se tienen los puntos de 


continuidad. 

El ordenador no puede, al menos en principio, intuir qué método a apli- 
car para estudiar la continuidad. La manera de averiguarlo es ir probando 
valores; cuando el ordenador, no el programa, dé error de operación ha- 
bremos encontrado un punto de discontinuidad. Lógicamente al interrum- 
pir el proceso, es necesario saber con qué valor ha ocurrido esto; para ello 
introducimos una instrucción PRINT antes de operar. Si se diese el caso 
de hallar uno, podríamos seguir buscando reanudando el proceso con lí- 
mite inferior, el valor de X siguiente al hallado. 

Este método rudimentario de buscar el crack del ordenador tiene en 
principio dos defectos: 


1.2  Obliga al operador a estar atento a la ejecución apuntando los pun- 
tos de discontinuidad. 

2.2 Puede ocurrir que aunque elijamos un incremento de X pequeño, 
el punto buscado puede estar entre dos valores de prueba, de tal 
forma que nos saltaríamos este punto y no sería detectado. 


Esta segunda dificultad puede presentarse aunque no se emplee el or- 
denador; por tanto, no es toda la culpa suya. 

De todas formas, a no ser que se analice previamente la función, es un 
método válido inicialmente. 
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=» >> 


80 LET X=X+(B-A)/L === 


=> TF X<B -FHEN GOTO 60 === 


Se ve cómo se puede elegir el intervalo a tratar, ya que los estudios de 
una función se hacen respecto a un intervalo concreto; sería imposible, 
quitando excepciones, un estudio exhaustivo de toda la curva. También se 
puede escoger el 0 incremento de X; a medida que el incremento sea in- 
ferior, mayor será la veracidad del estudio. 

Una vez estudiada la continuidad seguiremos con los demás conceptos. 
Tanto para máximo, mínimo, puntos de inflexión, etc., se necesita el con- 
cepto de derivada. 

La derivada de una función en un punto viene dada por: 


P(Xo)= lim f(Xo) —f(X+h) 
—P ——_—_—_—_——_—_——_—_— 
h 


Los límites, por el momento, no son asequibles al micro; para resolver 
este inconveniente se utilizará un método aproximativo. Partiendo de una 
H determinada, e ir decrementando este valor, se obtendrá una tabla de 
valores, de tal forma que llegará un valor de H en el que la diferencia de 
la derivada con la anterior hallada sea prácticamente nula. A partir de aquí 
se tomará como derivada el último hallado. Se utilizará como cota de la 
diferencia el valor 0,05, para no multiplicar en exceso los cálculos. 

Como siempre, para obtener la derivada de diferentes funciones, bas- 
tará modificar las instrucciones donde se aplique la función. Seguidamen- 
te se observará la estructura del proceso; la lista de las instrucciones da la 
subrutina que se empleará para el cálculo de la derivada en un punto. 


10 REM ESTUDIO DE LA DERIVADA EN. LUN PUNTO =— 
20 DEF FNA(X)=<Xx"2-1)- 
30 REM FPETICION DE NUMERO PARA CALCULAR SU DERIVADA. 
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Naturalmente, antes de probar un Xo habrá que estudiar su continui- 
dad, ya que ésta es condición necesaria para la existencia de la derivada. 

Conocido el proceso de la derivada pasaremos a estudiar la función. 
Empecemos viendo si una función es creciente o decreciente. 

Sea 


f0)=X?-1 
Estudiando el proceso en el intervalo [A, B] iremos comparando los dis- 
tintos resultados f(X), viendo cómo evolucionan. Se observará los distin- 


tos casos en procesos independientes, para no complicar el razonamiento 
(ver organigrama 20). 


La lista de instrucciones quedaría: 


e y ENTECIE! -. En 
S0-AF-EN=O-AFEN-—FRINT— ED EDITA IAN EL 


TS TIETPRES BEN IRINTIECIECTRI SMER ECRECIEN? 


270 T1E-CNSI- THEN PRINTES DECRECIENTE 
e zéá————_—_—_—_—___————— 


El uso de la variable de control es para definir si es estrictamente cre- 
ciente o decreciente. 

Ya se ha logrado avanzar otro paso para que el ordenador estudie una 
función. La dimensión de la matriz será (B—A)/L|, habrá que procurar 
que L sea divisor de (B—A) para no desperdiciar datos. 

Otros datos interesantes son los máximos y mínimos en un intervalo. Fí- 
jese en las gráficas siguientes: 


MAX 


Se observa que en un punto el valor de X es máximo o mínimo si la 
pendiente en ese punto de la tangente a la curva es nula. Para ello bastará 
que la derivada de la curva en ese punto sea 0. El sistema a seguir es es- 
tudiar la evolución de los valores de la derivada. Si cambia de signo querrá 
decir que existe un punto donde se anula la derivada. 

Ya se planteó en el capítulo 8 (la resolución de una ecuación) cómo se 
ve el punto de corte reduciendo el intervalo formado por un valor negati- 
vO y Otro positivo, o viceversa. 

Encontrando este punto habrá que averiguar si es máximo o mínimo; 
para esta labor, se observará la evolución del signo. 

Si la pendiente va de negativa a positiva, es un mínimo; en caso inver- 
so es un máximo. 
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Organigrama 20 
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Para observar la convexidad y concavidad de una curva hay que tener 
en cuenta los siguientes criterios: 


— Una función es convexa si f' es creciente. 
— Una función es cóncava si f' es decreciente. 


Si para ver los máximos y mínimos se estudió el signo de la derivada, 
para ver la convexidad y concavidad haremos un proceso similar al estu- 


dio de ver si una fuhción es creciente o decreciente al principio del tema. 
Ver organigrama 21. 


Organigrama 21 


Si no es cóncava ni convexa en el intervalo [A, B], quiere decir que exis- 
ten unos varios puntos en los que cambia; estos son los puntos de inflexión. 
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El programa para ver si es cóncava O convexa es: 


150 11€ ALCITi=aiiz4 THEN GATA 146 


1L£0 AF AtY)=AtY-1) THEN OIDO 190 


eel M ESTIMNIA 81 ES FANCAVE 
¿QO0 REM ESTUDIO 51 ES CONCAVA 


(ENAIYA)-ENALMO) Y) MA ME 
FNACALI=FNACXO?)7)/1xX1X0) 


ESA VET VAS INTIPERIA0O)/100 


I/O LES VA=INT(DERLODOO) /1000 


El lector puede intentar el cálculo de los puntos de inflexión. Para ello, 
deberá encontrar el punto donde cambia de creciente a decreciente, o vi- 


ceversa. ¡Inténtelo! 


LIZ 


AREA DEL RECINTO LIMITADO 
POR UNA CURVA. INTEGRALES 


I nos pidieran calcular el área limitada por una curva, a 
más de uno nos pondrían en un terrible aprieto. Rápida- 
mente intentaríamos hacer memoria de todas las reglas 
de integración, pero en la mayoría de los casos consegui- 
ríamos un terrible dolor de cabeza. 

Entonces la forma más fácil es, como anteriormente, 
pasárselo al ordenador. 

En principio, el ordenador «no sabe» integrar. Enseñé- 
mosle a hacerlo. 


dl 
] 
|] 
] 
u 
7] 
Wu 
a 


Si dividimos el intervalo en rectángulos iguales, la suma de estos rec- 
tángulos se aproximaría, más o menos, al área de la curva. Los rectángu- 
los exteriores darían un valor mayor que el de la curva y los interiores se 
quedarían «cortos». Si disminuimos el ancho de estos intervalos la dife- 
rencia se haría más pequeña. Repitamos este proceso varias veces consi- 
guiendo que el rectángulo sea prácticamente una línea, reduciendo la base 
de ellos; habremos conseguido una aproximación casi exacta al área de la 
curva. 
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En este momento, los rectángulos exteriores serían prácticamente igua- 
les a los interiores. 

Según este razonamiento, para «construir» nuestro programa debemos 
tener en cuenta: 


a) Cuáles son los límites de integración. 

b) Hemos de controlar razonadamente la disminución de las bases. 

Para ello dividiremos el intervalo completo en dos intervalos iguales, 
repitiendo progresivamente el proceso con intervalos que sean la mitad 
de los anteriores. De esta manera el primer cálculo sería con dos rectán- 
gulos, después 4, 6, 8 ..., así hasta conseguir que la longitud de los inter- 
valos sea prácticamente nula. 

c) La altura de estos rectángulos vendrá dada por el valor de la fun- 
ción en el extremo superior del intervalo. Nos estamos basando en los rec- 
tángulos exteriores para la resolución del problema; el razonamiento se- 
ría semejante con los interiores. Por tanto, nuestro programa valdrá para 
cualquier función, sólo bastará rectificar la instrucción donde se calcule 
la función. 


Para comprender mejor todo el proceso vamos a escoger una función 
sencilla, por ejemplo: f(X) =X? en el intervalo [a, b]. 

Ahora fijemonos en el organigrama y la lista de instrucciones siguien- 
te, inicialmente sólo lo que aparece en el interior de los recuadros 
sombreados. 

Según nuestro razonamiento, lo primero sería introducir los datos 
iniciales: 

INPUT “Introducir límites integración: ”; a, B. 


Entramos ahora en el proceso: repetitivo del cálculo del área. Hemos 
dicho que calcularíamos el área del recinto como suma de las áreas de los 
rectángulos circunscritos. Ahora bien, este proceso habrá que repetirlo su- 
cesivamente hasta conseguir una aproximación aceptable. 

Para obtener la suma en cada caso es necesario conseguir unos pará- 
metros básicos: 

116 SUMA=0 
120 J=3*2 


En la 116 inicializamos el valor de suma en la 120 y obtenemos el va- 
lor del intervalo, es decir, la longitud de cada una de las bases de los rec- 
tángulos. Observamos que J se duplica, luego I se reducirá a la mitad cada 
vez que iniciemos el proceso. 

Con la instrucción FOR, lo que hacemos es calcular el área de cada rec- 
tángulo e ir acumulando (SUMA) hasta obtener el valor del área total. Si 
nos fijamos, la altura (H) viene dada por el valor de la función en ese pun- 
to. Es aquí, en la 70, donde tendríamos que modificar para obtener la in- 
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tegral de cualquier otra función; la base, lógicamente, nos la da la longi- 
tud del intervalo (1) (ver organigrama 22). 


Introducir 
límites 
A, B 


ERROR 


I=INT(B-A)/J 


K=1 


¡ S X = -A+I+K 
AnremoR Ae 


SUMA + 
Mel 


NO 


SI 


sI 


Organigrama 22 


A E .- 
0 KEM CALCULO DEL AREA DEFINIDA FOR UNA CURVA 


20 INPUT “INTRODUCIR LIMITES DE LA INTEGRACION: ":48,EB 
O AA a AAC a 


ERROR=* "*C0 
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pus 


Una vez concluida la suma, imprimimos los resultados y volvemos a la 
115 a repetir el proceso. 

Se reiterará el proceso indefinidamente, ya que no hemos puesto nin- 
gún «control de fin». 

¿Cómo terminar? Podríamos tomar varias opciones, pero vamos a tra- 
tar dos. 

La primera de ellas (color verde) se basa en el hecho de que al aumen- 
tar J, disminuye I, llegando el momento que I fuese tan pequeño que la 
máquina lo considerase 0. En este momento, todas las áreas siguientes val- 
drían 0, luego el valor conseguido antes de este momento sería el buscado. 

Por ello, no tenemos más que preguntar si la SUMA=0. Terminando el 
programa en el caso de cumplirse. Podríamos haber hecho la pregunta des- 
pués de la XXX consultando el valor de I. 

La otra (color rojo) es poder nosotros elegir el final. Introducimos una 
cota de error y cuando la diferencia entre la suma y la suma anterior es 
menor que esta cota se termina el proceso. El porqué de esta opción se 
debe a que llega un momento que los rectángulos son tan estrechos que 
la diferencia área total es tan pequeña que no merece la pena seguir. 

Analicemos los resultados del programa: 


1 
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Calculamos ahora la integral por el método tradicional: 


1 1 
3 
fo [2] to 


Hemos comprobado que nuestro ordenador es capaz de hacer una in- 
tegral definida. Ahora bien, ¡cuidado!, ya que antes de hacerla habrá que - 
analizar esta integral. 

Debemos actuar como si la hiciéramos nosotros, es decir, ver los pun- 
tos de corte, continuidad... Nuestro ordenador es «listo», pero nunca po- 
drá hacer más de lo que le «enseñemos». 
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Se repasan los principales conceptos de la 
ciencia citada, desde un punto de vista 
eminentemente práctico y para su 
aplicación al ordenador personal. Se basa 
el texto en la presentación de pequeños 
programas (que usted podrá introducir en 
su ordenador personal); estos programas 
se analizan y se indican posibles 
modificaciones a introducir para variar 
los resultados que se deben obtener. Libro 
útil tanto para que el estudiante repase 


estos conceptos que está aprendiendo 
como para que el adulto recuerde algunos 
conocimientos olvidados que le interesan 
hoy. 


